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ABSTRACT 


Heat -transfer  measurements  were  made  for  a  single, 
water-heated  tube  in  a  pool  of  R-114  to  simulate  operating 
conditions  of  a  water  chiller.  Data  were  obtained  for  a 
smooth  copper  tube,  and  for  two  commercially  available 
tubes:  a  spirally  roped  copper-nickel  tube  with  a  porous 
coating;  and  a  copper  tube  with  a  structured  outer  surface 
and  a  multiple- start  helical  ridged  inner  surface. 
Measurements  were  made  for  refrigerant-oil  mixtures  at  oil 
concentrations  from  0  to  6  mass  percent  with  a  boiling  pool 
temperature  of  13. 8  °C.  Results  for  the  two  enhanced  tubes 
with  and  without  oil  are  compared  to  the  smooth  tube  data. 
Enhancement  factors  for  the  overall  heat- transfer  coeffi¬ 
cient  were  4.  0  and  3.  6  for  the  structured  surface  and 
porous- coated  tubes,  respectively,  in  pure  refrigerant  and 
at  a  water  velocity  of  2  m/s.  For  these  same  conditions  the 
enhancement  factors  A  for  the  outside  heatr-transfer  coeffi¬ 
cient  were  14.6  and  6.4  for  the  porous-coated  and  structured 
surface  tubes,  respectively. 
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Recent  developments  in  boiling  surfaces  have  shown 
considerable  enhancements  in  heat  transfer  performance.  The 
worldwide  literature  on  enhanced  heat  transfer  contains  over 
3000  published  technical  papers  and  reports  [Ref.  1:  p.  81]  . 
This  increased  interest  in  heat  transfer  augmentation  is  the 
result  of  incentives  for  energy  and  material  savings.  One 
effective  way  to  improve  the  heat  transfer  is  by  using 
passive  augmentation. 

Passive  augmentation  uses  fine-scale  alteration  of 
surfaces,  both  external  and  internal  [Ref.  1:  p.  82].  These 
surfaces  may  consist  of  either  an  applied  porous  coating  or 
fins  deformed  in  various  ways  to  provide  a  large  number  of 
reentrant  cavities  on  external  .  surfaces.  Methods  of 
internal  enhancement  include  insert  devices,  forged  fins,  or 
deformation  of  the  surface  (i.e. ,  spirally  roped  or  corru¬ 
gated  tube).  The  use  of  such  surfaces  can  lead  to  consider¬ 
able  reduction  in  the  size  and  weight  of  heat- transfer 
equipment.  The  reductions  may  lead  to  smaller  capital  costs 
or  operating  costs  or  both. 

Of  particular  interest,  in  naval  applications,  is  the 
reduction  in  the  size  of  water  chillers  in  refrigeration 
systems.  While  a  number  of  investigations  are  currently  in 
progress,  the  theoretical  treatment  of  the  boiling  perform¬ 
ance  of  various  tubes  is  almost  impossible  owing  to  the  very 
complicated  mechanisms  involved.  The  degree  of  difficulty 
increases  with  the  presence  of  oil  in  the  refrigerant 
liquid.  In  general,  refrigeration  systems  with  oil- 
lubricated,  hermetically  sealed  compressors  contain  small 
mass  percents  of  oil  in  the  refrigerant  liquid.  Therefore, 
reliable  data  covering  a  wide  variety  of  operating  condi¬ 
tions  for  various  refrigerants  and  different  boiling 
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surfaces  is  essential  to  further  develop  more  compact  and 
efficient  evaporators. 

Most  of  the  experiments  on  enhanced  boiling  surfaces 
reported  in  the  literature  have  used  electrically  heated 
tubes.  The  difficulties  associated  with  these  tubes,  espe¬ 
cially  during  the  instrumentation  stage,  raise  questions  as 
to  the  reliability  of  the  data  as  described  by  Wanniarchchi 
et  al  [Ref.  2:  p.  14]  .  Despite  the  precautions  taken  to 
minimize  contact  resistance,  it  may  still  be  present  between 
the  thermocouple  locations  and  the  outer  boiling  surface. 
Another  area  of  doubt,  described  by  Wanniarchchi  et  al. 
[Ref.  2:  p.  14] ,  is  the  nonuniform  heat  fluxes  generated  by 
commercially  available  electrical  heaters.  Additionally,  in 
order  to  use  an  electrical  heated  tube,  the  inside  surface 
of  the  tube  must  be  smooth,  or  any  internal  enhancement  must 
be  bored  out.  This  is  disadvantageous  as  designers  would  be 
interested  in  knowing  the  overall  heat- transfer  coefficient. 
Most  industrial  experiments  use  2-to  3-meter-long  tubes 
with  warm  water  flowing  inside  them.  While  the  information 
generated  from  these  experiments  closely  approximate  actual 
evaporators,  compared  to  electrical  heated  tubes,  the  large 
water  temperature  drop  from  inlet  to  exit  (up  to  5  K)  may 
make  it  difficult  to  study  the  details  of  the  boiling 
process. 

The  type  of  refrigerant  selected  will  also  influence  the 
design  of  evaporators.  R-114,  a  moderate-pressure  refrige¬ 
rant,  is  receiving  more  attention,  in  particular  for  naval 
applications.  Advantages  to  using  R-114  include:  (a)  it 
belongs  to  the  refrigerant  group  with  the  lowest  toxicity, 
(b)  it  is  very  stable  with  temperature,  and  (c)  it  has  a 
fairly  large  value  for  the  energy  transfer  per  unit  volume 
of  vapor  (i.e.  ,  Ev  =  h(g  /vg  )  [Ref.  3:  pp.  281-291].  The 
advantages  of  the  first  two  items  are  clear.  The  third  item 
is  advantageous  since  a  higher  value  for  Ev  means  lower 


pressure  drop  through  the  refrigerant  piping  for  a  specified 
heat  duty.  As  the  normal  boiling  point  of  the  type  of 
refrigerant  increases,  the  Ev  value  increases.  R-ll,  R-113, 
R-114,  and  R-22  have  Ev  values  in  increasing  order.  R-114 
may  be  preferable  to  R-ll  and  R-113  due  to  the  higher  Ev 
value  and  preferable  to  R-22  due  to  the  ability  to  use 
lighter  components  in  the  refrigeration  loop. 

Based  on  the  above  discussion,  the  major  objective  of 
the  present  investigation  was  to  study  the  boiling  perform¬ 
ance  of  two  commercially  available  tubes,  a  porous-coated 
tube  and  a  mechanically  structured  surface,  in  comparison 
with  a  smooth  tube  for  the  following  conditions:  (a)  tubes 
are  water-heated,  (b)  boiling  fluid  is  R-114  with  0,  1,  2, 

and  6  percent  by  mass  oil,  and  (c)  boiling  temperature  is 
13.8  °C. 


A.  NUCLEATE  BOILING 


When  the  temperature  of  a  boiling  surface  exceeds  the 
saturation  temperature  of  a  fluid  by  a  few  degrees,  nucleate 
boiling  occurs.  The  difference  between  the  boiling  surface 
and  saturation  temperature  is  the  amount  of  wall  superheat. 
Dougherty  [Ref.  4:  p.  175]  defined  pool-boiling  as  vaporiza¬ 
tion  occurring  under  the  following  conditions:  (1)  liquid 
depth  >>  the  bubble  diameter,  (2)  there  is  negligible  effect 
on  heat  transfer  due  to  low  or  no  externally  imposed 
velocity  of  the  fluid,  and  (3)  the  bubbles  move  away  from 
the  boiling  surface  due  to  a  force  field.  With  the  above 
conditions,  the  heat  flow  can  be  expressed  as: 


(2.1) 


A  description  of  the  boiling  mechanism  is  given  by 
Chongrungreong  and  Sauer  [Ref.  5:  p.  701].  A  thin  layer  of 
superheated  liquid  is  formed  adjacent  to  the  boiling 
surface.  Bubbles  nucleate  in  this  thin  layer  and  grow  from 
preferred  spots  on  the  boiling  surface.  It  is  assumed  that 
the  primary  resistance  to  heat  transfer  is  within  the  thin 
liquid  layer.  The  height  of  the  liquid  in  the  boiling 
container  is  not  a  primary  variable,  but  the  liquid  proper¬ 
ties  should  be  controlled.  However,  large  flooded  evapora¬ 
tors  may  experience  a  significant  "submergence  effect"  from 
the  liquid  head. 

3.  EXTERNAL  SURFACE  ENHANCEMENT 

For  smooth  heating  surfaces,  bubbles  nucleate  at  various 
scratches  and  cavities  on  the  surface.  Fujii  stated  that  the 
number  of  active  sites  increases  as  the  heat  flux  increases 
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[Ref.  6:  p.  48].  Webb  [Ref.  7:  p.  46]  stated  that  the 
ability  to  enhance  the  nucleate  boiling  coefficient  by 
applying  some  type  of  roughness  has  been  known  for  over  50 
years. 

Presently,  there  are  two  major  types  of  commercially 
developed  roughnesses:  (1)  porous  coating  and  (2)  various 
fins  and  surfaces  with  reentrant  cavities  [Ref.  8:  p.  24]. 
The  first  type  of  roughness  is  a  sintered  particle  coating, 
usually  copper  or  aluminum.  Figure  2. 1  shows  a  schematic  and 
a  photomicrograph  of  a  porous  coating.  The  small  reentrant 
cavities  are  interconnected  by  substrate  tunnels.  The 
surface  reduces  the  required  superheat  for  vapor  generation 
by  entrapping  a  high  density  of  relatively  large  vapor 
nuclei  in  the  cavities  contained  within  the  porous  coating 
[Ref.  8:  p.  24].  As  described  by  Webb  [Ref.  7:  p.  58], 
researchers  found  that  a  critical  pore  size  not  a  particle 
size  governs  the  number  of  nucleation  sites.  Large  pores  are 
required  for  liquids  with  high  surface  tension  and  high 
thermal  conductivity,  while  small  pores  work  best  for 
liquids  with  low  surface  tension  and  low  thermal  conduc¬ 
tivity  (i.e.  ,  refrigerants).  Czikk  and  O'Neill  [Ref.  9:  p. 
53]  developed  correlations  to  include  the  effects  of 
different  pore  sizes  of  porous-coatings.  They  concluded  that 
there  were  two  resistances  to  heat  transfer  for  these  coat¬ 
ings:  a  nucleation  superheat  related  to  bubble  diameter;  and 
a  conduction  superheat  controlled  by  the  liquid  film  sepa¬ 
rating  the  vapor  bubbles  from  the  metal  matrix.  Carnavos 
[Ref.  10:  pp.  106-108]  found  that  a  porous  coating  resulted 
in  a  700-800  percent  better  heat- transfer  coefficient  than  a 
smooth  surface  in  a  pool  of  R-ll.  Czikk  [Ref.  11:  p.  98] 
also  found  increased  performance  of  the  porous  -  coated  tubes 
used  in  his  20-ton  water- chiller  experiment. 

Structured  fins  or  surfaces  form  reentrant  cavities  of 
various  geometries.  Figure  2. 2  is  shows  an  example  of  a 
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-Reentrant  Cavity 
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.2.  Schematic  of  Manufactured  Reentrant  Cavity 
( iron  Ref.  3  ) . 
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manufactured  reentrant  cavity.  These  reentrant  cavities  act 
as  stable  nucleation  sites,  which  enhances  the  heat 
transfer.  The  nucleation  site,  to  remain  active,  is  depen¬ 
dent  on  the  mouth  diameter  falling  within  a  critical  range 
and  shape  with  a  maximum  reentrant  angle.  This  range  is  a 
function  of  the  fluid  properties.  The  cost  of  these 
surfaces,  presently,  are  not  significantly  higher  than  the 
cost  for  smooth  tubes  and  they  dramatically  improve  boiling 
performance  compared  to  the  smooth  tubes. 

The  mechanisms  by  which  the  reentrant  cavities  operate 
are  complicated.  Previous  research  concluded  that  the  combi¬ 
nation  of  bubble  evaporation,  thermal  boundary- layer  strip¬ 
ping  and  bubble  agitation  controlled  the  heat  transfer  from 
smooth  surfaces  [Ref.  12:  p.  192].  In  contrast,  different 
experiments  conducted  by  Arshad  [Ref.  12:  p.  192]  and  Arai 
[Ref.  13:  p.  37]  show  that  thin  film  evaporation  in  the 
reentrant  cavity  controlled  the  heat- transfer  mechanism. 
Bubbles  were  formed  by  vapor  exiting  cavities  as  the  liquid- 
vapor  interface  of  the  thin  film  contacted  the  cavity 
surface.  Surface  tension  holds  most  of  the  liquid  on  the 
cavity  walls.  Ayub  [Ref.  14:  p.  64]  showed  this  similar 
"thermosiphon  mechanism"  with  his  experiments  using  enhanced 
surfaces. 

Nakayama  [Ref.  15:  p.  37]  gives  a  fairly  detailed 
description  of  a  physical  model  undergoing  this  bubble 
growth  mechanism.  Figure  2. 3  shows  the  three  major  phases  of 
bubble  growth.  Phase  I  consists  of  a  pressure  buildup  in 
the  tunnel  by  evaporation  of  liquid  held  in  corners  and 
continues  until  the  meniscus  reaches  a  hemispherical  shape 
of  radius  r0  =  d0 / 2  (where  d0  is  the  mouth  diameter  of  the 
reentrant  cavity).  In  Phase  II,  the  meniscus  grows  faster  at 
some  pores  than  at  others.  The  Phase  I  pressure  buildup  is 
reduced  as  the  vapor  enters  the  growing  bubbles.  The 
meniscus  at  inactive  pores  does  not  grow  due  to  this  vapor 
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pressure  reduction.  Initially,  the  bubble  expands  under  the 
high  internal  pressure,  while  lacer  expansion  is  controlled 
by  the  receding  liquid  inertia  around  the  bubble.  Phase  III 
is  termed  the  liquid  intake  phase  with  liquid  flowing  into 
the  tunnel  through  inactive  pores.  This  flow  is  over  a  short 
interval  of  pressure  depression  which  occurs  as  the  pressure 
of  the  bubble  and  tunnel  is  lower  than  the  pool  pressure. 
The  bubble  leaves  the  pore  and  new  meniscus  formation  closes 
off  the  pore,  ending  this  phase  and  returning  the  cycle  to 
Phase  I. 


C.  EPFECTS  OF  OIL 

The  introduction  of  oil  into  the  pure  refrigerant,  in 
general,  causes  a  decrease  in  the  boiling  coefficient  as  the 
oil  concentration  increases.  The  effects  of  oil  on  the 
mixture  are  theoretically  complicated  as  most  empirical 
equations  relating  physical  properties  are  derived  form 
experimental  data  and  not  from  theoretical  equations. 
Various  empirical  equations  are  cited  by  Ghongrungreong  and 
Sauer  [Ref.  5:  pp.  703-705].  Yet,  Jensen  and  Jackman 
[Ref.  16:  p.  184]  showed  that  the  correlation  developed  by 
Chongrungreong  and  Sauer  appeared  to  overpredict  as  oil 
concentrations  increased  due  to  poor  prediction  of  the 
mixture  viscosity. 

Predictive  equations  were  developed  by  Jensen  and 
Jackman  [Ref.  16:  pp.  186-187]  for  density,  viscosity, 
surface  tension  and  specific  heat  for  pure  R-113,  pure  oil, 
and  R-113-oil  mixtures.  These  equations  are: 

,  .  C„ 


u 

.71 


(2.  3) 
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These  equations  are  based  on  a  physical  model  shown  in 
Figure  2.4  developed  by  Jensen  and  Jackman  [Ref.  16:  pp. 
187-188].  The  model  shows  a  bubble  growing  on  a  heated 
surface  as  the  refrigerant  evaporates  from  the  superheated 
liquid  phase  into  the  bubble  interior.  Due  to  the  less  vola¬ 
tile  nature  of  the  oil  in  the  mixture,  it  does  not  evaporate 
into  the  bubble,  but,  remains  at  the  liquid-vapor  interface. 
The  decrease  in  the  rate  of  bubble  growth  is  due  to  the 
decrease  of  oil  diffusion  into  the  liquid  mixture  and  the 
decrease  of  refrigerant  through  the  oil  layer  into  the 
bubble.  The  decrease  of  the  bubble  growth  rate  decreases 
the  heat- transfer  rate. 

D.  INTERNAL  ENHANCEMENT 

Internal  enhancement  increases  the  heat- transfer  area, 
creates  additional  turbulence,  and  produces  secondary  flows, 
all  of  which  contribute  to  the  increase  in  the  heat 
transfer.  The  laminar  sublayer  is  assumed  to  be  the  prin¬ 
cipal  resistance  to  heat  transfer  [Ref.  17:  p.  6].  The 
developed  secondary  flow  thins  the  sublayer  and  moves  it 
into  the  stream  where  turbulence  prevails  and  without  an 
increase  in  shear;  likewise,  heat  transfer  is  increased 
without  an  increase  in  friction  [Ref.  18:  p.  30].  Most 
studies  show  that  the  friction  factor  does  increase  with 
internal  enhancement  but  the  amount  is  strongly  dependent  on 
the  internal  enhancement  geometry.  The  increase  of  both  the 
friction  and  heat  transfer  is  dependent  on  fin  or  rib  pitch, 
groove  depths,  and  helix  angles  [Refs.  17,19:  pp.  6,  19-21]. 
Additionally,  these  studies  have  shown  that  there  are 


critical  geometrical  dimensions  that  give  the  maximum  heat 
transfer.  For  example,  short  fins  (<<  tube  diameter)  with 
rifling  will  raise  heat-transfer  coefficients  with  no  prob¬ 
lems  of  flow  stagnation.  Yet,  medium  fins  may  develop  flow 
stagnation  and  decrease  heat - trans f er  coefficients  [Ref.  19: 
pp.  22-23] .  As  such,  the  selection  of  internal  enhancement 
will  generally  improve  the  performance  of  boiling  surfaces. 


A.  OVERALL  APPARATUS 

A  schematic  of  the  experimental  apparatus  is  given  in 
Figure  3. 1.  The  photographs  in  Figure  3. 2  show  two 
different  views  of  the  apparatus.  Complete  details  on  the 
design  and  construction  of  the  refrigerant  and  oil  compo¬ 
nents  are  given  by  Karasbun  [Ref.  20:  pp.  24-32] .  A  discus¬ 
sion  of  the  modifications  for  the  water  heating  mode  is 
provided  in  this  paper.  The  major  components  of  the  appa¬ 
ratus  are:  two  Pyrex  glass  tees,  an  R-114  liquid  reservoir, 
a  water- ethylene  glycol  mixture  sump,  an  R-12  refrigeration 
system,  a  vacuum  pump,  a  water  supply  tank,  three  centri¬ 
fugal  pumps,  a  flow  meter,  three  heaters,  a  graduated  oil 
cylinder,  and  an  oil  reservoir.  The  R-114  boiling  and 
condensation  occurred  in  the  lower  and  upper  glass  tees, 
respectively.  The  R-114  vapor  was  condensed  by  the  water- 
ethylene  glycol  mixture  pumped  through  a  copper  condenser 
coil  located  in  the  upper  glass  tee.  A  1/2-Ton  R-12  refrig¬ 
eration  system  maintained  the  water-ethylene  glycol  mixture 
between  -18  and  -14  °C. 

The  major  difference  between  the  experimental  apparatus 
used  by  Karasabun  and  Reilly  [Refs.  20,8:  pp.  24-32,  30-35] , 
and  the  one  used  during  the  present  investigation  concerned 
the  use  of  water  heating  instead  of  electric  heating  as  a 
means  to  provide  heat  to  the  boiling  tube.  Filtered  tap 
water  supplied  from  a  storage  tank  (15),  was  pumped  by  two 
1/2-hp  motors  (13  and  14)  connected  in  series.  A  three-way 
valve  (V-16)  was  provided  on  the  downstream  side  of  the 
second  pump  to  obtain  two  flow  paths.  The  normal  flow  path 
for  data  runs  was  through:  the  metering  valve  (V-17),  the 
flowmeter  (6),  the  inlet  mixing  chamber  (16),  the  boiling 
tube,  the  outlet  mixing  chamber  (16),  the  chambers  provided 


with  heaters  (7  and  8),  and  returning  to  the  storage  tank 
through  valve  (V-20).  The  by-pass  flow  path  was  the  same  up 
to  the  three-way  valve  then  the  path  was  through  heater 
chamber  (8)  and  back  to  the  storage  tank  through  valve 
(V-20).  Valves  V-17  and  V-21  were  closed  in  order  to 
isolate  the  boiling  tube.  This  by-pass  option  was  essential 
for  preheating  the  water  to  higher  temperatures,  to  achieve 
higher  heat  flux  values  in  the  boiling  tube  compared  to 
values  obtained  from  room- temperature  water.  A  500-W  and  two 


1000-W  heaters  were  used  to  maintain  a  steady  water  inlet 
temperature  and  to  preheat  the  water.  Flexible  pressure  hose 
was  used  for  water  piping  throughout  the  apparatus.  The 
mixing  chambers  and  sections  of  hose  connecting  the  mixing 
chambers  to  the  boiling  tube  were  double  insulated  to  ensure 
accurate  temperature  measurements. 

Copper- constantan  thermocouples  were  used  to  measure 
thermal  emf ' s  for  R-114  liquid  and  vapor,  water  inlet  and 
outlet,  and  water-ethylene  glycol  mixture.  A  series- 
connected  thermopile  consisting  of  10  junctions  on  either 
end  was  used  to  measure  the  temperature  drop  of  the  water 
from  the  inlet  to  the  outlet  of  the  boiling  tube.  This  ther¬ 
mopile  was  calibrated  against  two  quartz  thermometer  probes 
and  agreement  was  found  to  be  better  than  ±  0. 02  K. 

B.  BOILING  TUBE  CONSTRUCTION 

Figure  3.3  shows  schematics  of  the  smooth  tube,  and 
two  enhanced  tubes  used.  The  tubes  tested  were: 

1.  a  smooth  copper  tube; 

2.  a  90:10  copper-nickel  corrugated  tube  (commercially 
referred  to  as  Korodense  tube)  with  an  external, 
sintered  porous  coating  (i.e.  High  Flux); 

3.  an  internally  and  externally  enhanced  tube  (alloy 
C12200)  produced  by  Wolverine  Division  of  U. 0. P. 

( coinmerically  referred  to  as  Turbo-B); 

4.  a  porous  -  coated  (High  Flux)  Korodense  tube  with  the 
porous  coating  machined  off; 

5.  a  Turbo-B  tube  with  the  external  enhancement  machined 
off. 
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Schematic  of  Boiling  Tubes:  fa)  Smooth  Tube 
(^bj  High  Flux  Couteu  Tuoe,  and  (c)  Turbo-3 
Tube. 


These  tubes  will  be  referred  to  as  High  Flux,  Turbo-B,  modi¬ 
fied  High  Flux,  and  modified  Turbo-B  tubes,  respectively, 
with  the  exception  of  tubes  2  and  4.  Tubes  2  and  4  will  be 
referred  to  as  Korodense  tubes  with  respect  to  discussions 
of  only  the  internal  enhancement.  All  of  the  tubes  had  an 
active  (i.e.  ,  heated)  section  of  304.8  mm  in  length.  The 
smooth,  High  Flux,  and  modified  High  Flux  tubes  had  an  outer 
diameter  of  15. 9  mm  and  an  inside  diameter  of  12. 7  mm.  In 
the  case  of  the  High  Flux  tube,  these  diameters  represent 
the  values  in  the  smooth  portion  of  the  tube  (i.e. ,  before 
performing  the  corrugation  process).  The  Turbo-B  tube  had 
an  external  structured  surface  and  a  multiple-start  helical 
ridged  inner  surface.  This  tube  had  an  outer  diameter  of 
16. 9  mm  at  the  base  of  the  external  structured  surface  and  a 
minimum  diameter  of  14.  5  mm  at  the  internal  ridge  tips.  The 
modified  Turbo-B  tube  had  identical  diameters  as  the  Turbo-B 
tube.  The  63.5-mm-long  inactive  sections  on  either  end  was 
insulated  by  Teflon  sleeves  of  1. 6  mm  wall  thickness  located 
inside  the  boiling  tubes,  as  shown  in  Figure  3.3.  The 
requirement  for  using  these  two  modified  tubes  is  explained 
further  in  Chapter  IV  (EXPERIMENTAL  PROCEDURES). 

C.  COMPUTER- CONTROLLED  DATA  ACQUISITION  AND  REDUCTION 

Hewlett-Packard  equipment  was  used  for  data  acquisi¬ 
tion,  control,  and  analysis,  as  well  as  for  storing  of 
data.  An  HP-3497A  acquisition  unit  was  used  to  read  the 
thermocouple  and  thermopile  outputs.  Table  1  gives  the 
channel  designations  used. 

Information  entered  by  keyboard  to  an  HP-9826A  computer 
unit  prompted  and  controlled  the  data  acquisition  unit. 
Data  were  analyzed  and  stored  with  the  same  computer  unit.  A 
step-by-step  description  of  the  data- reduction  procedure  is 
given  in  Appendix  A,  along  with  a  printout  of  the  program 
( DRP4 )  used. 


TABLE  1 

CHANNEL  DESIGNATIONS 


Channel 

0-1 

2 

3 

4 

5 
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Function 

R-114  liquid  thermal  emf's  T(l) 
and  T(2) 

R-114  vapor  thermal  emf_T(3) 


water- ethylene  g 
thermal  emf  T(4; 


lycol  mixture 


water 

water  _ 

thermal  emf 
water  inlet 


inlet  thermal  emf  T(5 
outlet  thermal  emf  T( 


difference 
and  outlet 


h 


across 
T(  7) 
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A.  FLOWMETER  CALIBRATION 


This  section  gives  a  brief  description  of  the  flowmeter 
calibration;  a  more  detailed  explaination  is  given  in 
Appendix  B.  A  Fischer  Porter  flowmeter  was  used  to  indicate 
percent  water  flow  through  the  boiling  tube.  Prior  to 
conducting  any  boiling  runs,  it  was  necessary  to  develop  a 
correlation  relating  the  flow  percent  to  the  mass  flow  rate 
and  water  velocity.  Additionally,  any  temperature  effect  on 
viscosity  would  be  incorporated  into  this  correlation  as  a 
correction  factor.  Calibration  runs  were  conducted  over  a 
temperature  range  of  19  °C  to  38  °C  and  varying  flowmeter 
settings.  The  correlation  (equation  (B. 1))  was  determined 
from  these  results.  The  correction  factor  for  the  viscosity 
temperature  dependency  was  determined  to  be  1.  0  with  an 
accuracy  of  ±  0.02  as  explained  in  Appendix  B.  In  fact,  it 
was  not  possible  to  find  a  systematic  trend  for  this  correc¬ 
tion  factor  with  the  water  inlet  temperature.  The  stated 
accuracy  mainly  consists  of  the  uncertainties  involved  with 
the  visual  reading  of  the  flowmeter  settings. 

B.  BOILING  TUBE  AND  APPARATUS  PREPARATION 

The  external  surface  of  the  boiling  tube  was  cleaned 
with  a  Nitol  (2%  nitric  acid  and  98%  ethyl  alcohol)  solution 
and  both  the  external  and  internal  surfaces  were  rinsed  with 
acetone  and  were  air  dried  prior  to  installation  in  the 
apparatus.  A  vacuum  test  at  an  absolute  pressure  of  about  50 
mm  Hg  was  conducted  on  the  refrigerant  side  of  the  apparatus 
after  the  boiling  tube  was  installed.  If  no  signs  of  leaks 
were  evident  after  two  hours,  the  system  was  pressurized  to 
a  gage  pressure  of  about  250  mm  Hg  with  R-114  vapor  and  the 
apparatus  was  checked  for  leaks  with  an  automatic  Halogen 


Leak  Detector.  After  fixing  any  leaks  that  may  have  been 
present,  R-114  liquid  was  transferred  from  the  reservoir  to 
a  pre-determined  level  in  the  boiling  tee.  This  level  gave  a 
free  surface  40  mm  above  the  centerline  of  the  boiling  tube. 
At  the  set  level,  the  mass  of  the  R-114  liquid  was  computed 
to  be  2.  48  kg.  A  desired  saturation  temperature  of  4  °C  and 
the  resultant  system  pressure  were  obtained  by  turning  on 
the  R-12  refrigeration  system  and  varying  the  flow  of  cold 
(approximately  -16  °C)  water- ethylene  glycol  mixture  through 
the  condenser  coil  with  valve  V-9.  At  this  saturation 
temperature  and  pressure,  each  boiling  tube  developed 
nucleation  sites.  The  saturation  temperature  of  4  °C  was 
maintained  for  at  least  30  minutes  prior  to  initiating 
nucleate  boiling  on  the  complete  active  length  of  the 
boiling  tube. 

The  initiation  of  nucleate  boiling  could  be  done  with 
cold  water  at  19  °C  or  with  warm  water  at  25  °C  through 
valve  V-17  or  with  steam  from  a  steam  generation  system 
through  valve  V-22.  Early  experimental  data  showed  slight 
differences  (up  to  10  percent)  in  the  computed  outside  heat- 
transfer  coefficient  depending  on  the  highest  heat  flux  at 
which  the  boiling  was  initiated  on  each  boiling  tube.  As 
will  be  discussed  in  Chapter  V  (RESULTS  AND  DISCUSSION),  the 
steam  intiation  method  provided  the  highest  initial  wall 
superheat  and  ensured  the  be^t  possible  repeatability  of  the 
data.  This  intiation  method  was  used  for  all  of  the  runs 
unless  otherwise  specified. 

Steam  initiatior  was  conducted  by  introducing  steam 
through  valve  V-22  and  discharging  steam  drain  valve  V-18 
for  about  one  minute.  After  this  period,  the  three-way  valve 
V-16  was  shifted  from  a  neutral  position  to  the  normal  flow 
position.  Valves  V-22  and  V-18  were  closed  and  valve  V-21 
was  opened  simultaneously  as  the  three-way  valve  was  shifted 
to  the  normal  position.  Valve  V-17  was  opened  wide  prior  to 


shifting  of  the  three-way  valve.  As  these  valves  were  shift¬ 
ed,  warm  water  at  25  °C  was  pumped  through  the  boiling  tube 
at  the  maximum  water  velocity,  which  in  turn  gave  the  maxi¬ 
mum  heat  flux  for  the  given  water  inlet  temperature. 

C.  BOILING  DATA  RUNS 

Once  the  nucleate  boiling  initiation  was  completed,  the 
saturation  temperature  was  raised  to  13.  8  °C  by  adjusting 
valve  V-9.  This  valve  was  also  used  to  maintain  this  satura¬ 
tion  temperature  within  ±0. 05  K  throughout  a  boiling  run. 
The  water  inlet  temperature  was  maintained  to  within  ±0. 08  K 
by  using  the  heaters  or  introducing  filtered  cold  water, 
respectively.  The  combination  of  heaters  and  cold  water  was 
used  to  maintain  the  desired  water  inlet  temperature  during 
each  run.  Two  types  of  boiling  runs  were  conducted  at  each 
oil  mass  percent. 

The  first  type  of  run  was  conducted  after  holding  the 
desired  saturation  temperature,  water  inlet  temperature  and 
initial  flowmeter  setting  at  steady-state  conditions  for 
about  10  minutes.  After  this  period,  the  flowmeter  setting 
was  decreased  with  two  data  sets  being  taken  at  each  flow 
setting.  The  two  data  sets  were  to  show  repeatability  at 
each  setting.  The  period  between  pairs  of  data  sets  was 
about  5  to  8  minutes,  with  steady  state  for  each  flowmeter 
setting  maintained  for  at  least  2  minutes.  The  second  type 
of  boiling  run  consisted  of  maintaining  the  previously 
mentioned  initial  steady-state  conditions  at  a  selected  high 
water  inlet  temperature.  The  run  was  commenced  by  intro¬ 
ducing  cold  water  into  the  storage  tank  which  lowered  the 
water  inlet  temperature.  While  holding  the  water  velocity 
constant,  six  to  seven  data  points  were  taken,  on  a  one  time 
pass,  at  about  every  0. 3  K  inlet  temperature  decrease.  When 
the  inlet  temperature  reached  a  selected  low  temperature, 
the  steady-state  conditions  of  the  low  inlet  temperature. 


water  velocity,  and  saturation  temperature  of  13.8  °C  were 
maintained  for  about  10  minutes.  Then,  the  water  inlet 
temperature  was  increased  using  the  heaters,  while  holding 
the  water  velocity  constant.  Again,  data  points  were  taken 
on  a  one  time  pass,  in  about  0.3  K  intervals,  until  the 
inlet  temperature  reached  the  previously  used  high  inlet 
temperature. 

Two  complete  boiling  runs  of  the  first  type  mentioned 
above  were  performed  at  0%  and  2%  oil  concentrations  to 
further  demonstrate  repeatability.  Table  2  shows  the  percent 
difference  in  wall  superheat  (i.e.  ,  AT)  and  heat  flux  for 
each  boiling  tube  at  these  two  oil  percents.  The  higher 
percentages  for  the  High  Flux  and  Turbo-B  tubes  for  the  wall 
superheats  are  simply  because  of  the  low  wall  superheat 
values  of  these  two  tubes  when  compared  to  the  values  for 
the  smooth  tube.  At  these  low  superheat  values,  a  small 
difference  in  the  wall  superheat  leads  to  a  larger 
percentage  difference  than  that  which  occurs  at  higher  wall 
superheats.  The  variation  of  the  water  inlet  temperature  was 
not  considered  a  factor  in  this  percent  difference  as  the 
variation  was  less  than  ±  0.02  K  per  minute  during  these 
constant-inlet-temperature ,  decreasing-water-velocity  runs. 

TABLE  2 

BOILING  RUN  REPEATABILITY 


Tube  0il%  A  T 


q 


Smooth  0 
High  Flux  0 
Turbo-B  0 
Smooth  2 
High  Flux  2 
Turbo-B  2 


4 

±30% 

±30$ 

±10$ 

±10% 

±30% 


j2l? 
±  2% 


Note:  initiations  performed  with  steam 


Additionally,  data  runs  were  conducted  to  determine  the 
temperature  increase  across  the  boiling  tube  due  to  the 
pressure  drop  between  the  locations  of  the  thermopile 


probes.  The  data  collected  showed  that  this  difference  was 
less  than  the  ±0.02  K  accuracy  of  the  thermopile  and  was 
considered  negligible. 

Following  the  boiling  runs  in  pure  R-114,  oil  was 
introduced  to  the  boiling  tee  through  valve  V-l  from  the 
graduated  cylinder  (5)  shown  in  Figure  3.1.  After  the 
required  volume  of  oil  for  a  given  oil  mass  percent  was 
added,  water  at  the  maximum  velocity  was  pumped  through  the 
boiling  tube  to  promote  vigorous  boiling.  This  boiling 
ensured  good  mixing  of  the  R-114  and  oil  mixture.  Either 
another  boiling  run  was  conducted  or  the  system  was  shut 
down  to  prepare  for  another  steam  initiation. 

D.  DATA  ACQUISITION  AND  REDUCTION 

The  automatic  data  acquisition  system  was  prompted  to 
record  the  required  thermal  emf ' s ,  by  keyboard  inputs  to  the 
computer.  The  data  were  immediately  processed  and  printed 
out  on  a  hard-copy  printer.  A  step-by-step  procedure  of  data 
processing  and  reprocessing  is  given  in  Appendix  A.  Also, 
a  printout  of  the  data  reduction  program,  DRP4,  is  included 
in  this  appendix.  The  heat  flux,  overall  heat-transfer  coef¬ 
ficient,  outside  boiling  heat- transfer  coef f icient ,  and  wall 
superheat  values  were  computed  based  on  the  outside  area 
expressed  by  the  diameter  if  the  enhancement  of  the  outer 
surface  was  removed.  Further,  in  order  to  account  for  the 
heat  conduction  through  the  inactive  end  sections,  a  correc¬ 
tion  was  made  for  heat  flux  using  an  iterative  computation 
procedure  based  on  natural  convection  as  discussed  by 
Karasbun  [Ref.  20:  pp.  54-56]. 

When  a  boiling  run  was  completed,  the  data  were  repro¬ 
cessed  by  computer  to  obtain  the  inside  and  outside  heat- 
transfer  coefficients  for  further  calculations.  A 
Sieder-Tate- type  constant  (Cj  )  for  the  inside  coefficient 
for  each  of  the  externally  smooth  tubes  was  calculated  by 
using  a  modified  Wilson  plot  |^Ref.  22]  .  For  the  smooth  tube, 
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modified  High  Flux,  and  modified  Turbo-B  tubes,  the  outside 
coefficient  was  calculated  simultaneously  using  the  cor¬ 
relation  (equation:  (C.  3)  and  (C.4))  developed  by  Rohsenow 
[Ref.  23:  p.  969]  ,  and  the  modified  Wilson  plot.  The 
Rohsenow  correlation  is  based  on  externally  smooth  tubes, 
thus  the  requirement  for  the  modified  High  Flux  and  Turbo-B 
tubes.  A  further  explaination  of  the  modified  Wilson  plot 
is  given  in  Appendix  C.  The  C|  values  found  for  the  exter¬ 
nally  smooth  tubes  were  used  for  the  corresponding  exter¬ 
nally  enhanced  tubes.  Knowing  the  inside  coefficient,  the 
outside  coefficient  for  the  boiling  tubes  were  computed  by 
subtracting  the  inside  and  wall  resistances  from  the  meas¬ 
ured  overall  thermal  resistance  (equation:  (C. 7)). 


A.  INSIDE  HEAT-TRANSFER  COEFFICIENT 


Results  were  obtained  for  the  smooth,  High  Flux,  and 
Turbo-B  boiling  tubes  at  a  boiling  temperature  of  13. 8  °C 
and  at  oil  concentrations  of  0,  1,  2,  and  6  mass  percent. 

The  heat  flux  versus  wall  superheat  for  these  three  tubes  is 
shown  in  Figure  5. 1.  All  three  tubes  are  in  the  nucleate 
boiling  range  with  decreasing  heat  flux.  It  can  be  seen  that 
the  High  Flux  and  Turbo-B  tubes  outperform  the  smooth  tube 
throughout  the  tested  heat  flux  range.  Additionally,  the 
Rohsenow  correlation  for  a  smooth  tube  at  a  boiling  tempera¬ 
ture  of  13.8  ,°C  with  the  experimentally  determined  (from  the 
modified  Wilson  plot)  CSf  coefficient  of  0.0060  is  shown  in 
this  figure. 

The  Rohsenow  correlation,  as  used  in  the  modified  Wilson 
plot,  has  an  exponent  of  r  =  1/3.  Yet,  this  correlation  is 
very  sensitive  to  the  r  value.  A  change  in  this  exponent  of 

three  percent  yielded  a  change  in  the  CSf  only  of  0.  016 

percent  but  gave  a  change  of  17  percent  in  the 

Sieder-Tate- type  constant  inferred  from  the  modified  Wilson 
plot.  For  the  purpose  of  data  analysis,  the  value  or  r  =  1/3 
was  used,  while  knowing  that  an  additional  uncertainty  was 
being  introduced  by  tins  exponent. 

For  each  externally  smooth  tube,  the  Sieder-Tate- type 
cons",  ant  was  computed  using  the  modified  Wilson  plot  based 
•  n  '  "o  1  ;  different  runs  in  pure  R-114.  The  computed 
vi.  : s  t  r  "he  smooth,  Korodense,  and  Turbo-B  tubes  were 
:  1.  ').0o6  t  0.005,  and  0.077  ±  0.003,  respec- 

"  F-  r  ".lie  smooth  tube,  the  experimentally  found 

dioder  Tate  type  constant  is  33  percent  greater  known  than 
the  well  known  value  of  0.027  for  long  internally  smooth 


appears  to  be  the  result  of  the  entrance  effects  of  the 
shorter  experimental  tube  and  the  uncertainty  introduced  by 
the  exponent  r  of  the  Rohsenow  correlation  as  discussed 
above.  The  minimum  entrance  length  for  a  fully  developed 
pipe  flow  is  given,  as  a  rule  of  thumb,  by  Incropera  and 
DeWitt  [Ref.  24:  p.  406]  as 
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for  Reynolds  numbers  greater  than  or  equal  to  10000.  Also, 
as  the  pipe  roughness  increases  the  minimum  entrance  length 
decreases. 

Withers  has  developed  Stanton  number  correlations  for 
tubes  with  single-helix  and  multiple  helix  ridging 
[Refs.  25,26:  pp.  52-56,44].  These  correlations  are  listed 
below.  For  single-helix  ridging, 


St 


./r-0'5  (i 

l)0*53 

7.22 

(!)  M 

rj 

a.  127 

♦  Y 

l  * 

—• 

(5.2) 


where  r  and  m  are  determined  by  tube  dimensions,  and 
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For  multiple-helix-  ridging, 
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and  m  again  are  functions  of  tube  dimensions, 
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A  comparison  of  Stanton  numbers  was  made  for  the  modified 
Korodense  and  modified  Turbo-3  tube3  using  the  following 
equations  for  these  two  tubes: 
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The  results  for  these  two  modified  tubes  were  compared  to  a 
copper  Korodense  (i.  e.,  single-helix  ridging)  and  a 
Turbo-Chill  (  i.  e.,  multiple-helix  ridging)  tubes,  respec¬ 
tively.  The  latter  two  tubes  tested  by  Withers  were  similar 
in  radial  dimensions  but  were  1.  5  meter  long.  The  predicted 
St  numbers  using  equations  (5.1)  -(5.8)  for  the  Withers’ 
Korodense  and  Turbo-Chill  tubes  were  0.00177  and  0.00273, 
respectively.  The  experimental  St  numbers  using  equations  ( 5.9 ) 
and  (5.10)  for  the  modified  Korodense  and  modified  Turbo-3 
tubes  were  0.00242  and  0.00286,  respectively.  The  majority 
of  the  difference  between  the  predicted  and  experimental  St 
numbers  may  be  attributable  to  the  entrance  effects  and  the 
uncertainity  introduced  by  the  Rohsenow  correlation. 

3.  LIGHT  OF?  EFFECTS 

As  referred  to  in  Chapter  IV,  Figure  5.2  shows  the  light 
off  effects  for  the  smooth.  High  Flux  (i. e.,  porous- coated 
Korodense),  and  Turbo-B  tubes.  These  effects  are  most  prob¬ 
ably  due  to  the  requirement  of  greater  heat  flux  to  initiate 
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the  nucleation  sites  than  the  heat  flux  required  to  maintain 
the  sitfjs  once  they  are  activated.  As  the  boiling  tube  sits 
in  the  1-114  liquid,  prior  to  initiation,  reentrant  cavities 
or  the  nucleation  sites  flood  with  liquid.  The  initial 
formation  of  the  vapor  bubbles  and  thus  the  initiation  of 
the  nucleation  site  requires  the  input  of  a  greater  heat 
flux  or  higher  wall  superheat.  Figure  5. 3  shows  the  boiling 
tube  initiation  with  steam  for  the  smooth,  High  Flux,  and 
Turbo-B  tubes  [Refs.  21,7:  p.  478,  62].  As  shown  in  Figure 
5.3,  the  High  Flux  tube  is  more  sensitive  to  light  off 
effects.  This  sensitivity  is  mostly  due  to  the  lower  wall 
superheats  at  which  nucleation  occurs.  A  change  of  . 5  of  a 
degree  is  more  pronounced  at  these  lower  wall  superheats 
than  a  similar  change  at  the  wall  superheats  for  the  Turbo-B 
and  smooth  tubes. 

C.  SMOOTH  TUBE 

Figure  5.4  shows  the  performance  of  the  smooth  tube  at 
0,  1,  2,  and  6  mass  percent  oil  concentration  at  13.8  °C. 
Reilly's  [Ref.  8:  p.  65]  data  for  the  electrically  heated 
smooth  tube  at  a  boiling  temperature  of  -2.2  °C  in  pure 
R-114  liquid  are  also  included  for  comparison.  The  differ¬ 
ences  between  the  experimental  results  shown  and  Reilly's 
data  are  probably  due  to  physical  characteristics  introduced 
in  the  electrically  heated  tube  when  soldering  the  copper 
sleeves  (inside  of  which  the  electric  heaters  were  fitted) 
on  the  interior  of  the  tube.  As  discussed  by  Reilly  [Ref.  8: 
pp.  57-61] ,  the  contact  resistance  at  the  interface  between 
the  sleeve  and  the  inner  surface  of  the  tube  was  minimized 
by  tinning.  However,  the  tinning  process  may  not  have  been 
100  percent  successful,  thereby  introducing  unacceptable 
uncertainties  into  the  experimental  measurements.  Further 
Reilly  reported  considerable  variations  (up  to  3  K  at  a  heat 
flux  of  98  kW/m2 )  in  the  measured  wall  temperatures.  Reilly 
attributed  this  observation  to  an  axially  non-uniform  heat 
flux  generated  by  the  cartridge  heater. 


As  shown  in  Figure  5.4,  the  wall  superheat  increased  for 
1  percent  oil  concentration  and  decreased  slightly  for  2  and 
6  percent  oil.  This  decrease  in  the  wall  superheat  for  the 
higher  oil  percentage  is  possibly  due  to  the  enhanced  bubble 
formation  experienced  by  the  tube  because  of  the  foaming 
action  of  the  R-114  and  oil  mixture  decreasing  resistance  to 
heat  transfer  as  explained  by  Chongrungreong  [Ref.  5:  p. 
701]  ,  Nobukatsu  [Ref.  27:  p.  60]  ,  and  Chaddock  [Ref.  21:  p. 
474]  . 

D.  HIGH  FLUX  TUBE 

The  performance  of  the  High  Flux  tube  in  terms  of  heat 
flux  versus  wall  superheat  is  given  in  Figure  5. 5. 
Additionally,  this  figure  shows  Reilly's  data  [Ref.  8:  p. 
70]  ,  for  the  electrically  heated  High  Flux  tube  in  pure 
R-114  at  a  6. 7 °C  boiling  temperature.  The  presence  of  oil 
increased  the  wall  superheat,  at  a  heat  flux  of  40  kW/m2 ,  by 
a  factor  of  1. 22  for  1  percent  oil  and  a  factor  of  2.  12  for 
6  percent  oil.  The  number  of  nucleation  sites  reduced 
considerably  and  the  foaming  action  increased  significantly 
as  the  oil  mass  percent  increased.  These  actions  are  shown 
in  Figure  5.  6.  With  6  percent  oil  concentration,  it  was 
visually  observed  that  the  nucleation  sites  were  active  only 
in  the  corrugations,  as  shown  in  Figure  5.6  (c).  This 
congregation  of  sites  is  due,  in  part,  to  the  sparser  porous 
coating  on  the  high  points  of  the  tube  and  the  decreased 
wall  thickness  in  the  corrugations.  As  a  result  of  the 
corrugation  process  during  tube  manufacture,  the  wall  thick¬ 
ness  at  the  corrugation  is  smaller  than  the  rest  of  the 
tube,  which  implies  lower  wall  resistance.  During  the 
sintering  process,  the  coating  appears  to  have  concentrated 
in  these  corrugations  giving  a  thicker  coating  of  copper 
particles  and  increased  number  of  nucleation  sites  compared 
to  the  rest  of  the  tube  surface. 


The  increase  of  the  wall  superheat,  when  oil  is  present, 
is  possibly  the  result  of  increased  surface  tension  and 
viscosity.  The  increasing  surface  tension  and  viscosity  lead 
to  the  formation  of  an  oil-rich  layer  next  to  the  boiling 
surface.  This  oil-rich  layer  has  an  insulating  and  inhib¬ 
iting  of  bubble  growth  effect  within  the  cavities.  The 
higher  liquid-vapor  surface  tension  and  surface  tension  in 
the  oil-rich  layer  requires  more  energy  for  the  formation 
and  growth  of  the  bubble.  Additionally,  diffusion  processes 
affect  the  bubble  formation.  Oil  concentrations  near  the 
boiling  surface  are  high  compared  to  the  concentration  in 
the  bulk  liquid  and  diffusion  occurs  from  the  cavity  into 
the  bulk  liquid.  This  diffusion  continues  until  an  equilib¬ 
rium  point  is  reached.  If  oil-rich  layers  are  developing  on 
the  boiling  surfaces,  the  diffusion  from  the  cavities  to  the 
bulk  liquid  slows  down  aiding  the  further  build-up  of  the 
oil-rich  layer.  The  vapor  bubble  forming  in  the  cavity  also 
has  an  oil  vapor  content  which  is  undergoing  its  own  diffu¬ 
sion  process.  As  the  refrigerant  evaporates  through  the 
oil-rich  layer  into  the  interior  of  the  bubble,  it  must 
overcome  large  diffusion  resistances  set  up  by  the  oil  vapor 
within  the  bubble.  The  increased  surface  tension,  viscosity, 
diffusion  resistance  and  the  decreasing  bubble  formation 
rate  contribute  to  the  decrease  in  the  heat- transfer  coeffi¬ 
cient  and  the  increase  in  the  wall  superheat  at  a  given  heat 
flux.  [Refs.  5,21,28,29,30:  p.  702,  477-478,59,372,82] 

Also,  boiling  runs  were  conducted  holding  the  water  vel¬ 
ocity  constant  and  allowing  the  water  inlet  temperature  to 
decrease  and  then  to  increase.  The  decreasing  and  increasing 
heat  flux  with  the  respective  change  in  the  water  inlet 
temperature  versus  wall  superheat  is  shown  in  Figure  5.7. 
The  motion  of  the  hysteresis  ( i. e. ,  clockwise  or  counter¬ 
clockwise)  depended  on  the  temperature  starting  point.  If 
the  temperature  was  decreased  and  then  increased,  the  motion 


was  clockwise.  When  the  temperature  was  increased  and  then 
decreased,  the  motion  was  counterclockwise.  This  directional 
change  is  related  to  the  amount  of  heat  flux  required  to 
maintain  active  nucleation  sites.  In  the  first  case  of  tem¬ 
perature  decreasing  then  increasing,  less  heat  flux  is  need¬ 
ed  to  maintain  the  same  amount  of  sites,  but,  at  the  same 
time,  heat  flux  decrease  as  the  water  heat  capacaity  is 
decreasing.  In  the  second  case,  with  the  water  inlet  temper¬ 
ature  increasing  then  decreasing,  the  water  heat  capacity 
increases  so  the  heat  flux  increases  raising  the  wall 
superheat  and  the  heat  -  transfer  coefficient. 

E.  TURBO -B  TUBE 

Figure  5. 8  shows  a  similar  relationship  as  Figure  5. 5 
but  for  the  Turbo-B  tube.  The  increase  in  wall  superheat, 
for  a  heat  flux  of  40  kW/m2 ,  ranges  from  a  factor  of  1.04 
for  1  percent  oil  to  a  factor  of  1. 37  for  6  percent  oil.  The 
Turbo-B  tube  remained  active  over  the  full  active  length  of 
the  tube,  as  shown  in  Figure  5.9.  The  apparent  decrease'of 
nucleation  sites,  as  shown  by  the  High  Flux  tube,  did  not 
occur  with  the  Turbo-B  tube.  As  discussed  pertaining  to  the 
High  Flux  tube,  the  increase  in  the  wall  superheat  with 
increasing  oil  concentration  is  probably  due  to  the 
increased  surface  tension  and  viscosity  of  the  mixture 
within  the  cavities.  Also,  Boiling  runs  for  constant  water 
velocity  and  varying  water  temperature  were  performed  on  the 
Turbo-B  tube  and  are  also  shown  in  Figure  5. 10.  The  same 
hysteresis  motions  observed  with  the  High  Flux  tube  are 
observed  with  the  Turbo-B  tube. 

F.  COMPARISON  OF  BOILING  HEAT-TRANSFER  COEFFICIENTS 

The  ratios  shown  in  Figure  5. 11  are  h/h  for  the 
Turbo-B  and  High  Flux  tubes  and  h/hs  for  the  smooth  tube. 
All  of  the  ratios  were  taken  at  a  heat  flux  of  40  kW/m2.  The 
High  Flux  tube  shows  the  most  significant  effect  due  to  oil 
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concentration  which  is  due  to  the  reasons  previously 
discussed.  The  h/hs  ratio  for  the  High  Flux  tube  decreased 
by  65  percent.  For  the  Turbo-B  tube  the  rate  of  degradation 
of  the  boiling  heat-transf er  coefficient  is  less  pronounced. 
The  Turbo-B  tube  shows  a  decrease  of  24  percent  at  6  mass 
percent  oil.  The  larger  cavity  size  of  the  Turbo-B  tube, 
compared  to  the  High  Flux  tube,  results  in  smaller  oil 
concentrations  in  the  oil  rich  layer.  This  observation  may 
explain  the  smaller  degradation  experienced  by  Turbo-B  tube. 

G.  OVERALL  HEAT-TRANSFER  CHARACTERISTICS 

The  overall  heat- transfer  coefficient  versus  decreasing 
water  velocity  for  oil  mass  percents  of  0,  2,  and  6  for  each 
boiling  tube  is  shown  in  Figure  5.12.  For  pure  R-114  with  a 
water  velocity  of  2  m/s,  the  Turbo-B  tube  outperforms  the 
High  Flux  tube  by  a  factor  of  1. 13  and  the  smooth  tube  by  a 
factor  of  4.  While  the  High  Flux  tube  outperformed  the 
Turbo-B  tube  based  on  the  outside  heat- transfer  coefficient, 
the  reverse  is  true  for  the  overall  heat- transfer  coeffi¬ 
cient.  This  relationship  in  U0  is  due  in  part  to  the 
increased  internal  enhancement  of  the  Turbo-B  tube. 
Comparison  of  the  Sieder-Tate- type  constants  for  the  Turbo-B 
and  High  Flux  tubes  (see  p.  34)  shows  that  the  Turbo-B  tube 
has  a  17  percent  greater  inside  coefficient  than  the  High 
Flux  tube.  Both  of  the  enhanced  tubes  dramatically  outper¬ 
formed  the  smooth  tube.  It  should  be  noted  that  the  increase 
in  pressure  drop  due  to  internal  enhancement  (compared  to  a 
smooth- interior  case)  was  not  considered  here.  The  pressure 
drop  would  be  an  important  factor  in  finally  selecting  the 
most  economical  tube  type. 
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A.  CONCLUSIONS 


Considering  the  data  gathered  in  this  investigation  for 
boiling  of  R-114  and  oil  mixtures  at  13.8  °C  boiling  temper¬ 


ature,  the  following  conclusions  are  reached: 


1. 

2. 


3. 


4. 


Based  on  the  overall  heat- transfer  coefficient,  the 
Turbo-B  and  High  Flux  tubes  outperformed  the  smooth 
tube  by  a  factor  greater  than  3,  and  the  Turbo-B  tube 
outperformed  the  High  Flux  tube  by  a  factor  of  1.  13. 


At  a  heat  flux  of  40  kW/m2  with  zero  mass  percent  oil, 
the  outside  coefficients  of  the  High  Flux  and  Turbo-fi 
tubes  were  factors  of  14.6  and  6.4  times,  respec¬ 
tively,  compared  to  the  smooth  tube.  These  factors 
decreased  to  7.0  and  4.9  for  the  High  Flux  and  Turbo-B 
tubes,  respectively,  with  6  mass  percent  oil. 


The  Turbo-B  and  High  Flux  tubes  showed  Sieder-Tate- 
type  constants  of  0.077  and  0.066,  respectively,  com¬ 
pared  to  the  Sieder-Tate  constant  of  0. 027  for  long 
smooth  tubes.  As  noted  in  the  discussion,  these  values 
may  be  up  to  10  percent  lower  for  internally  enhanced 
long  tubes  than  the  values  for  the  short  tubes. 


While  the  High  Flux  tube  outperformed  the  Turbo-B 
tube,  based  on  the  outside  heat- transfer  coefficient, 
the  High  Flux  tube  is  more  susceptible  to  the  presence 
of  oil. 


B.  RECOMMENDATIONS 

Based  on  the  results  obtained  from  this  investigation, 
the  following  recommendations  are  made: 

1.  The  investigation  of  water-heated  tubes  should  be 
upgraded  to  use  long  tubes.  This  would  decrease  the 
entrance  effects  that  short  tubes  experience  and  it 
would  come  closer  to  duplicating  actual  operating 
conditions.  Also,  the  pressure  drop  effect  previously 
mentioned  in  Chapter  V  (RESULTS  AND  DISCUSSION}  could 
be  investigated  in  the  analysis  of  the  long  tubes. 

2.  Studies  of  the  boiling  heat- transfer  coefficients 
should  be  expanded  to  tube  bundles.  The  interaction  of 
tubes  within  the  tube  bundle  is  a  major  factor  in  the 
analysis  of  heat- transfer  data.  The  boiling  heat- 
transfer  coefficient  may  be  strongly  influenced  by  any 
foaming  action  within  the  bundle. 

3.  Further  investigations  should  be  conducted  on  internal 
tube  enhancements  due  to  the  major  thermal  resistance 
being  the  internal  resistance. 

4.  Modifications  to  the  condensing  side  of  the  present 
apparatus  would  allow  for  studies  of  short  water- 
heated  tubes  at  lower  boiling  temperatures.  The  major 


modification  would  involve  only  the  increase  of  the 
condenser  drain  piping  to  handle  increases  in  the  flow 
rate  of  the  refrigerant  liquid. 
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DATA  REDUCTION  PROGRAM 


The  data  reduction  program,  DRP4,  used  for  this  investi¬ 
gation  is  listed  below.  A  brief  description  of  the  major 
sections  is  given  in  Table  3,  which  is  followed  by  an  indi¬ 
vidual  line  listing  of  DRP4. 


TABLE  3 

DRP4  MAJOR  SECTIONS 


Line  Number 


Description 


10-80 

80-3690 


3700-4165 


4170- 

4225- 

4940- 

6180- 

6415- 

6495- 


4220 

4935 

6175 

6410 

6490 

7280 


7285-7390 

7395-9640 

9645-9725 

9730-9755 

9760-10200 


10205-10360 

10365-10475 


Selection  options  (i.e.  process 
data ,  plot , etc.  ) 

Sub  Program  Main 
-boiling  tube  dimensions 
-select  desired  operating 
conditions  (i.e.  initiation 
mode,  flow  rate,  Tsat,  etc.) 
-convert  emf  s  to  temperatures 
-compute  contact  resistance 
(electrically  heated  mode  only) 
-compute  various  water  properties 
-various  heat  transfer 
calculations( Q,  LMTD,  U  ,  etc.) 
-compute  various  Freon  114 
properties 

-compute  natural  convective  heat- 
transfer  coefficient  for 
unenhanced  ends( iterative 
procedure ) 

-compute  heat  loss  for  unenhanced 
ends 

-compute  actual  heat  flux  and 
boiling  coefficient 
-record  and  store  data 
Functions  for  curve  fits  of 
various  Freon  114  properties 
Function  for  polynomial  curve  fit 
Sub  Program  Poly 
-calculate  nth  order  polynomial 
Sub  Program  Plin 
-plot  linear  graphs 
Sub  Program  Stats 
-calculate  standard  deviations 
Sub  Program  Coeff 
-plot  cross-plot  files 
Sub  Program  Wilson 
-calculate  modified  Wilson 
plot  values  (see  Appx.  C) 
Functions  for  curve  fits  of 
various  water  properties 
Sub  Program  Plot 
-plot  non-linear  graphs 
Sub  Program  Symb 
-print  various  symbols  on 
plotted  graphs 
Sub  Program  Purge 
-purge  unwanted  files 
Sub  Program  Tdcn 
-calculate  water  temperature 
rise  corrections  due  to 
pressure  drop  between 
thermopile  locations 
Sub  Program  U  plot 
-print  U  files 
Sub  Program  Select 
-select  various  sub 
program  options 


-»c. 

ail 

^5 

50 

50 

55 

71 

75 

30 

55 

‘30 

35 

Too 

~nS  < 
*  ''5 

’  -o 

t  IS 

I  °0 

•25 
T  30 
’25 

1  40 
’  42 

150 
'55 
1  60 
1  63 
170 
1  75 
•30 
'35 
"?0 


:o5 
1 ' 5 


RE'/ISZO: 

. :  20/ 
PRINTE5  IS  ' 

OAI _ Seiect 

INPUT  "NAMT 


i  r 

r  :  .  .1 / 


^NQTHE^  jPtI3N  (  t  *  y  ,  0«N) .  Ise  i 


Ir  Isei-'  'HEN  OO'C  40 
3EE° 

3EED 

PRINTER  IS  ' 

PRINT“"DATA  COLLECTION/REPROCESSING  COMPLETED” 

END 

SUB  •iin 

COM  /  lao/  IdD 

COM  /Co/  C<7>.Icai 

COM  /'rii  ;/  DO  .D  i  ,  Do  .L  . Lu. <cu 

DIM  Ent  <  12>  .  T  t  12>  .  D 1  a*  3  >  .D2a<3)  .  D  ia*  3 )  .Doa<^?>  .  !_a(9 )  .  Lua<  3  >  .  Xcuat  9  > 
4) : '5  : 

DATA  •).  ' 003609 1 .25707 .34353 . -7S7345 .3095 . 78025535 .31 
DATA  -3247 486589 . b . 37633E * ’ 1 , -2 .661 32E-1 3 . 3 . 340 7BE*1 4 
BEAD  CO 

’  DATA  "Smoo  th”  .  ”H  t qn  r  I llx ”  .  "Thernoexe  L -E”  .  "Thermoexel -HE" 

DATA  Smooth, Hiqn  r  lux  .  TurDo-3  .H  ign  Eiux  Mod.Turoo-3  Mod 
READ  7n5(*) 

°9 INTER  IS  701 
3EEP 

IF  IdD-4  'HEM  PRINTER  IS  1 
IF  Iao-4  THEM  GOTO  1230 

INPUT  "ENTER  MONTH.  DATE  AND  rIME  ( MM : DD :HH :MM : SS> ” . DateS 
OUTPUT  703 : "'D" ; Dates 
OUTPUT  709 : "TD" 

ENTER  709: DateS 
PRINT 

PRINT  "  Month,  date  and  time  DateS 

3RIN' 

PQI‘lT  USING  "*0X. ""NOTE:  Program  name  :  DRP4 . 

3EE3 

INPUT  "ENTER  DISK  NUMBER". Dn 

PRINT  USING  "1SX.""Di;k  numoer  -  "".ZZ":Dn 

3EE3 

INPUT  "ENTER  INPUT  MODE  <0“305*»6. 1  •FILE)’*.  In 
3E£p 


INPUT 


HEATING  MODE  <0=ELEC;  1 -HATER)' 


[hm 


:co 

3EE= 

c 

INPT  "ENTE 

R  "HER 

MOCOUPLE  Tvpc 

<n*NEW. i-uLD>". 

leal 

‘  a  c, 

IF  ln*0  THE 

'■‘5 

3EE? 

:so 

INPU'  "S'h'E 

A  NAME  -?R  THE  RAW 

data  PILE". DO 

i  1  eS 

'55 

PRIM'  USING 

“1  6X  . 

•"‘Ne-ii  ‘lie  name  :  ""  .  1  4«"  ;  02 

i  le'3 

6  0 

3 1  :e ' » 20 

ORE  A"  3D AT 

D2^;i 

ie  5 . S i 2  9 ; 

-  C  1 

ASS  I  ON  »F  .  . 

9  U_  'J 

d:_*  :  ;*•* 

'?  1  ' 

i  ^1  ;;<m v  r  *;  r 

jnt:l 

'Inin  <MOUM 

I'.'C 

5*  ■  .  . eS~"D 

1  Jf*4M  V  •• 

•  >> 

IREA’E  3CHT 

V  •*. 

!  a*  5  .  J  ;  :e  : 

>  c 

393  I  ON  'P  .  . 

'0 

:<  ’9) 
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IN PL'  "EN'ER  lUNBER  OF  DEFECTIVE  ’C3  '  0  *DEF4UL  ' > " . Id t c 
i7  i  a  t  c  *  o  then 

_  2  t  C  !  *  ') 

j  * ,-  ■■ .  n 

PRINT  is  IMG  defect:  iv»  ’Cs  *xist . 

EMC  :7 

Ir  Ictc”:  *HE,N 
3EED 

IMP'JT  "ENTER  DEFECTIVE  TC  LOCATION”  .Ldtc ' 

PRINT  USING  "1bX."*’TC  is  detective  at  location  .  D“ :  Late  ' 

_jtc2-0 
END  IF 

IF  Idtc”2  THEN 
3E£° 

INPUT  "ENTER  'DEFECTIVE  rC  LOCATIONS". Ldtc: .LdtcE 

PRINT  USING  "T6X."”TC  are  defective  at  vocations  "”  .  D  .  '<X  .  3"  ;  Ldtc 1  .  .J  tcE 
END  IF 


AGO 

IF  Idtc>2  ’HEN 

405 

3EE° 

A 1  0 

PRINTER  IS  1 

A  1  5 

3EE? 

A2Q 

PRINT  "INVALID  ENTRY' 

AOS 

PRINTER  IS  "01 

430 

GOTO  310 

A35 

END  IF 

A  A!) 

END  IF 

445 

OUTPUT  iP;  le: :Ldtc1  . 

Ldtc. 

A50 

?  I«i”1  oorion 

A55 

ELSE 

AGO 

3EE? 

AG5 

INPUT  "GIVE  THE  NAME 

OF 

A70 

PRINT  USING  "16X.—0 

Id  f 

a  75 

ASSIGN  »Fiie2  TO  32_ 

•lie 

ABO 

ENTER  «F i Ie2 : Nr an 

-*B5 

ENTER  3F : .  e2 : Do i aS 

3RINT  USING  " 1 6X . ""T 

his 

A?5 

ENTER  »F : le2 : Ld tc ‘  _ 

3 1  cE 

500 

IF  Ldtc ! >0  OR  _atc2>0  TH 

'  4A" ; Do  i  o5 


PRINT  USING  "'6X,""rhonocoucies  iiere  cere  cti’/e  at  iocar  ions  :  ""  .  2  <  3D  .  JX )  ’ 
:  ■  .v-atc2 
END  IF 

ENTER  4F. L*2: Itt 
END  IF 

'  IF  I, '!»••)  AND  Ihn* !  'HEN  '535 
3EE= 

INP'jt  "NANT  '0  CREATE  A  PLOT  FILE ’  ( 0-N  Isiot 

IF  Iptot”;  THEN 

3EED 

INPUT  "OI’-'E  NAME  PAR  PLOT  pIl£”.P 
CREATE  3DA  ‘  p_-'  .  .eS.- 
ASS  IGN  :>A  .  D  t  'I  J_-  ;  .-5 

end  :r 

I-  I  nm” 1  'HEN 

rcr: 

INP1,’ '  "nan T  'C  IREA 'E  7ILE^  •*'»".  I 


o-  «> 


Eft 


:  rt 

'Tit:; 

.o  - 1  _  E 

'iH'4E"  .  .  -  *?s 

7’* 

**J  *  .  .  t?  _■  .*  r  ;  - 

3 

*5  '  '  ' 

END  Ir 

n  '  5 

CEE3 

620 

INPUT  '•'.-.AM* 

"G  CREATE 

Re  FILE'’  <  0-N  .  ’  *  '  .  I 

525 

’  *-  r  o  •»  '  J 

■530 

3EE3 

-  I MP'J T  "ENTER 

Fe  -ILF 

NANE"  .  .Ref  :  .  e  3 

541) 

CREATE  3DA'  : 

3  t  i  7 '3  . 1 

o 

5-i5 

ASSIGN  -Ref; 

le  TO  Re f 

:  1*3 

650 

END  IF 

•D3D 

END  :r 

660 

DR INTER  IS 

DO  0 

IF  In-0  THEN 

572 

3EE’ 

575 

PR INT  OS  IMG  ' 

4 X .  Seiect  taoe  numoer 

680 

IF  Ihm-0  ’HEN 

555 

pRINr  OS IMG 

''5X .  '"'0 

smooth  4  men  Ret 

530 

PRINT  USING 

"5X  . . 

Smooth  4  men  Cu  <Press/S 

535 

PRINT  JSING  ' 

"GX .  ‘"‘2 

Sot  t  soieer  4  men  >,a 

700 

PRINT  USING 

"EX. ■••'2 

Sott  Soiaer  4  men  HIGH  r 

705 

PRINT  USING 

"5X .  ‘”*4 

Hieiand  Hard  3  men"'*" 

7’  n 

PRINT  USING 

"GX .  '"'5 

HIGH  FLUX  3  men . 

7  ’  5 

PRINT  USING 

"GX .  "”5 

Gc’HA-K  19  r  ms/ in 

720 

PRINT  USING 

"GX.""7 

GEHA-<  26  “ms/in . 

725 

PRINT  USING 

"GX . . 3 

GE^tA-r  13  Eme/;n . 

730 

PRINT  JSING 

-EX. “"2 

GEAA-T  2S  F  i ns/ : n . 

755 

ELSE 

740 

PRINT  JSING 

•’GX  .""0 

Smooth  tuce  ’ 

745 

PRINT  USING 

"GX. . 

Hign  - .ax 

750 

PRINT  JSING 

"GX .  '"'2 

| aroo'3 

755 

PRINT  USING 

"GX . "“5 

H  ign  r  lax  Nod 

750 

PRINT  USING 

"SX.— 4 

Turno-3  Noa 

’53 

END  :r 

INP'JT  Itt 

-»  -»  ^ 

f)U-P’Jr  iF.ie 

•  ;  l  :  t 

7  ■?  3 

END  IF 

",'i5 

DR~'\*ER  12  " 

1 

’■?o 

~~  7  -  *  <  ’  *1  'H 

EN  30  IN  r 

JSING  "’GX  .  ",,Taoe  Nunoer: 

"°5 

IF  I :  t > 6  ’’HE 

N  POINT 

JSING  " '  6X  .  "iioe  NanDer: 

Z  Z  ."i  :  oa  ,  iioe  ;  .l/u  .  ..  ;.  * 

:.-  :w  "-eu  =ri.nt  using  "*5x.— Tod©  :yP@:  ••••.  isa";  Tns<  itt> 


3  '  'j 

3  -  - 

:.NP»jr  * 
FFrD 

•ENTER  OUTPUT 

VERSION 

<  0 -luNG . 1 -SHQR T . 2 -NONE ) " , lov 

d-  n 

3'c  1 

INPU  ”  ' 

'TELE  I'  O-.IO 

.  *  wpp  .  2 

3^LI'S^hap>/2)",  1 1  <3  v 

•)  .  .  ■  *  1 
325 

>  Di‘4E‘ic; 

•  r  _  r  -  *  ^ 

::ms  -s  'e:~ee 
17  -e^te:  mose 

*UBES 

1  ;  .■*  | 
D^'J 

'  3  *  * 0  i  djne r e r  it  t  n«?  r^o-roup  i  2 

3os i z  ions 

345 

la-A  . 

25. .01 ‘1 

125.  .,1123540,  .  0 1  2  44  A  .  .0123540.  .01  00355 

JC  j 

ic: 

-•  u 

3r.;r  j 
I '  =* .  i 

J*:nr?65.  .  J 1  *  E  F 

*  i‘  *  • 

<  *  -  *  > 

. . J1 157 

jo* rut s . te  diameter  };  .inennancea  arcs 

3ft 'A  .  0 ' 5575  .  .3*5375.  .3*5375.  .3*582*.  .0*5375.  .3’ 532*.  .3:  270.  .3:1 
READ  Doa<  *  > 

-•'.jriot .i  s»  annanced  sur*ace 

DATA  >0;s. . tots.  .  :015. . 1016. .2022. .2022. .2032. .2032. .2032. .2032 
READ  _a>  * > 

_u"'-eno  .-.  of  inennanced  sur-acs  at  the  ends 

DATA  .025*. .025m. .025*. .025*. .0752. . 0752 .. 0762 .. 0762 , .0762. .0762 
READ  Laa(*' 

<cu*76ernai  Conduct  iv i  tv  of  tuoe 

DATA  258 . 3** .3** . *5 . 2** . *5 . 3** , 9** . 398 . 298 

READ  <caa(*) 

IF  Ihm*  ‘  THEM 


Data  statements  for  jater  beating  node 


005’ 

on 

DATA  1.0 ’5375 

on 

READ  D2a<  * ) 

020 

DATA  1.0127.0 

025 

READ  D ia<  *  > 

030 

DATA  0 . J ! 5675 

035 

READ  Doa‘ *  > 

0*0 

DA 7  A  0.20*8.3 

0*5 

RE*0  _a<  * 1 

550 

DATA  1 . 025* . 3 

055 

READ  _  j.a<  *  > 

060 

DArA  393. *5. 2' 

065 

READ  *.cua<  *  > 

071 

EMD  IF 

;t  -c 

D2*D2a< 1 1 1 ) 

no 

D  1  *D  i  a'  1 1 1 ) 

135 

Do  *Qoa<  1 1 1  > 

no 

_  -  _  a.  t  I  tt  > 

095 

La*aaa< 1 1 1  > 

no 

<ca*'d ua<  1 1 1  > 

•  nc 

X  n  *  .3 

’ -o 

”r*  .  3 

•  •  c 

IF  ltt*0  ’HEN 

*  20 

IF  Itt>9  'HEN 

•  25 

*-* I *<  Do  2-Di 

•  'h  n 

O.P  T .Da 

’55 

3EE6 

*  *^0 

IMP’J T  ••■'JJ9E  I! 

IF  It  m-1  "HE 

*50 

IF  It  m-2  'HE 

*  c  r 

H  *■  i  _r^  3  ’►*  C* 

*  *D  <J 

IMP1:'  *  "'2,',P  'V! 

n  itv-o  'hem 

■  — 

I  -  ~  _  '*E  1 

."in-1  HEM  -niNT  ■  J  S  TNG  bX  .  ""  *"ixDe  Initiate: 


EM  3DIMT  JSIMG 


■  UDe  .nitiate: 
'ice  Initiate: 


*ot  water 

5  team . 

Co  id  Mater' 


IMP'."  '"E'*P'VE-_  MODE:  <  0  -  *  -.CMS  7  . ;  -DEC : :  *  r  -  DEC  .  >  -CONS  T :  2*7- INC .V-C 


;as i ng/ Mo ns  tan  r 


i:  ;'  v  j: 

-22 

Ir  I:t*2  IP  I-:-4 

N  ' . 3 . 

2  •' ' 

'NP':'  "EN  'ER  ’.I  i  3Er  : 

»H*  .  222  .  -“f  -  .  O' 

*  Q 

PRI.V  USING  "5a.  ""2. 

e<:er-  a:-? 

220 

3RI.NT  .SING  '"5X."" 

lonsranc 

C 

rND  Ir 

::u 

Ir  Inn*'  AMD  In*’  HMD 

.  j  i  1  3  !  "EN 

IF  1 1 1  *  )  "WEN  I .  *  .  022 

.2  JO 

IF  1 1 1 ■ ’  ]R 

N  2. * . J53 

'u  5 

IF  I:t*2  ]R  1 1 1  *  a  'HE 

N  C.- .382 

250 

ASSIGN  *F;;eL  TO  * 

255 

2AI _ A  1 1  son  <  C  r  .  2  ;  > 

2S0 

ASSIGN  »F ;  1  e2  '0  22  ■* 

:  !e5 

255 

ENTER  *F ;  ie2  *. Mrun  .  Do  i. 

dS ,_dtc1 ,Ldtc2 

2") 

END  IF 

Ns 'ID  *  0 

230 

IF  I do  *  a  them  Ihm*1 

2^5 

IF  Inm*1  THEN  Nsud*o 

230 

N  t  c  *  5 

^  a  ’ 

200 

Sx-0 

305 

Sy  *  0 

310 

3xs*0 

3  ’  5 

S  X  V  *  0 

Reoeat :  ' 

1225  IF  Im-0  THEM 
1  220  3 1 id*  2 . 22 

'325  Iao*2 

•2oi)  ;]N  <EV  3.15  RECOVER  '220 

•R45  PRINTER  *S  ' 

•25.:  ?9lNT~'jSI>tG  "ox  .  '‘"SELECT  DPT:CN . 

'255  PR  INI  USING  "SX . ""0* TAKE  DATS . 

1  250  IF  ISm-0  'HEN  PRINT  USING  "SX. -GE*  HEAT  CL’JX . 

' 365  IF  Ihm-'  'HEN  PRINT  J$ING  "SX  . •  $£*  WATER  FLOW  RAT' 

1270  2N T  USING  "5X ,  ”'2*SE'  ";at . 

■  2 75  PRINT  USING  "AX .""NOTE :  KEV  1  -  ESCAPE . 

' 230  3EE? 

•  235  INPUT  Ido 

•  290  IF  Ido>2  "HEN  Ioo*2 
'225  IF  I  do  *  0  than  2 ‘ j5 

'  41)!)' 

'405'  .OOP  '0  GE"  hE.AT  rL'JX  2R  FLO  WMF  "ER  SE'TING 

•  a  *  2  IF  I  do  * ‘  "HEN 
*-i  ’  5  IF  1 7i n ■  •)  r E ' I 

'^2')  5ur^jr  7Q9:MAR  AF*2  ^L*j  V°5*‘ 

'  ^20  INP'JT  "ENTER  DESIRED 

3RINT  USING  "AX  .  "”2ES IPFS  God  AC T; jAf 


iCD 


-i<3CD»‘-'0  1  <  3  I'OE"1-  > 

Ir  -iB3  <  SqcD-ijooB  i  ">£ r r  " 
IF  -taaD>0caD  '-iEM 
3EED  4000..: 

3EE3  4000 . .2 


•520 

3EE3  4000.  .2 

'525 

;loc 

'530 

3EED  250.. 2 

1535 

3EE?  250.. 2 

■  5^0 

3EEP  250.  .2 

1545 

END  IF 

•550 

PRINT  USING  ‘ 

!  555 

WAT  2 

'560 

G0T3  t445 

'565 

ELSE 

15  70 

3EE3 

•575 

PRINT  USING  ' 

’  5 30 

Err *500 

1535 

<iflr  2 

•530 

GOTO  ’445 

t  59c; 

END  IF 

1800 

ELSE 

1805 

3  EE? 

1810 

INPUT  "ENTER 

•6-5 

GOTO  1350 

1820 

END  IF 

•825 

END  IF 

•620' 

'  525  1 

!  LOOP  'U  3E’ 

’640 

IF  Ido-2  'HE' 

’  tHS 

IF  Ikat*'  7  HI 

i  S50 

3EE° 

’  555 

INPUT  "ENTER 

"560 

!  PRINT  USING  1 

’  SoS 

Ikat* 1 

"67  0 

3  i  a  ’  *0 

1675 

3ia2*0 

’  630 

Nn*  " 

•635 

Nrs-Nn  NQD  " 

1630 

Nn  •Nn**' 1 

’635 

IF  N  ra  *  1  ”WE' 

'  700 

PRINT  USING 

1  oat 

■705 

END  IF 

•  *  .'i 

IF  IFm  =  'l  *'~E' 

•7-5 

IF  I hm *  *  'HE 

•  720 

~3R  >•  "U  5 

«  ?  ■'VC 

IF  Ihm-0  61-ID 

’  730 

Sujn*0 

■7  35 

OUTPUT  703 : " 

’  7  4f! 

-3R  *0 

■  7 

ENTER  '99 :£. 

j  an  3  : 

‘  '55 

NE:<*  -  : 

■'b>, 

-  .  :  a  * lim  /  2  0 

.  I  *"  \  /S1/  »  ■ 

•  ■*  *•  c 

■=>  ;-*£ '  * 

•  -o.- 

I-  >: 

DTsat  ATsat  Rate 


;:»'jr  ’03:  "sr  ar3  rlh  vrs" 


»X! 

ffi 

•i| 


w 


•  -  ;c 

IF  .  -  J  -s.  1  sumo- 

*  7 'Vi 

3 ^  -it 0  ”  •  n  ^  * 

‘  735 

IF  I "S'  "HEN  '.jus-') 

'  Sip 

nex”  : 

"305 

IF  16m-'  ~HEM 

i  a:  o 

0UTDIJ”  709  :  "OR  AF29 

•  3  1  5 

OUTPUT  709 : "PS  OPT 

'  320 

Sum  *  0 

"325 

FOR  "0  20 

1330 

ENTER  709  :E 

1335 

Sum-Sum -E 

18^0 

NEXT  <lc 

1345 

Emt ( 7  > -4gS( Sum/20  > 

’850 

To  1 ie-Emt ( 7 ) /3 . 96E~ 

1855 

END  IF 

1360 

A  t  i  d  -  <  T I  d  1 -Tld2)*.5 

1865 

IF  ABS < fl t ia-D 1 id > >  . 

1370 

IF  Atid>Dtlo  THEN 

1375 

BEE?  4000  .  .2 

1330 

9EE?  4000.. 2 

’885 

BEE3  4000 . .2 

1890 

ELSE 

1895 

BEEP  250  .  .2 

1300 

SEE3  250.  .2 

1905 

BEE3  250.. 2 

toil) 

END  IF 

1915 

Err  1 -At la-Oldl 

1920 

0  i  d  1  -  8 1  i  d 

«  925 

Err2*  7v'-01d2 

’930 

Gld2-w 

•925 

IF  Ihm-0  'HEN  3RINT 

1940 

IF  Ihm-’  AND  lap*0 

UiDP  .  . 

iniet . Toi ie 

•  345. 

IF  Ihm-1  AND  Ido-4 

si  .  *1 

d2  .  Tv  .  'sump .‘in 

•950 

HAIr  2 

1  955 

GOTO  ’635 

f  '350 

ELSE 

•965 

IF  ABS<  A t id-D t id  >  > . 

'970 

IF  9tid>0tid  'HEN 

"375 

BEEP  3000  .  .2 

•330 

BEE3  3000  .  .2 

’995 

ELSE 

•990 

BEEP  300 . .2 

•  395 

3EE?  3 00.. 2 

2000 

END  IF 

2005 

Err ’ -At ;a-]la1 

20  ‘  0 

0  I  d  1  »  A  t ;  d 

20:5 

Err2*'v-']id2 

2020 

Gld2-'v 

2025 

IF  Ihm-0  "HEN  PRINT 

2030- 

IF  Ihm- 1  'HEN  PRINT 

.  0  iiiTio  .  :  n  i*? ,  ,di.j, 

20^0 

GGT0  '6.35 

7  ft  h  g 

ELSE 

:o5o 

3EED 

■>  <•»  c  - 

E  r  r  *  ^  v -  0  1  - 1 

-  s  0 

lidi-j-.d 

2057 

!  . T ld2 . Tv , T samp 
" ;  0 1  id  .  T  lal  .  7  \aZ  . ' 

,3(t13D.DD.2X>";Dti 


I  ,  I  1 Oc  >  i  V  , 

'  X  >" ;Dt  io , 


muwiP'PW 


257  1  ir  i *■  t ' "  -no  i-.n-o  '-E'i 

'630  ”«•<  '!"'))*■(•;))  .'2 

2635  END  IF 

2630  IF  Ihm-0  THEN  27 ’0 
2595  IsujnD* 7 <  1  2'Nsud  > 

2700  7iniet-*(13~Nsuo> 

2705  Tout  -  *  < '  4-N'sud  ) 

27 ’0  IF  Ihm-0  THEM 
27*5  Amo-085 < I r  > 

2720  No  1  t*OBS(Nr ) *25 
2725  Q-Volt-Amp 
2730  END  IF 

2735  IF  Itt-0  AND  Ihm-0  THEN 
2740  Kcu-FNKcu<  T.ja) 

2745  ELSE 

2750  Kcu-<cua( Itt) 

2755  END  IF 
2750  ! 

27S5!  FOURIER  CONDUCTION  EQUATION  HITH  CONTACT  RESISTANCE  NEGLECTED 
2770  IF  Ihm-0  THEN  Tui-7.aa-<>'-0G<  D2/D1  >  /  ( 2*PI*<cu-L  > 

2775  IF  I  lav-0  THEN  Tsat-Ttaa 

2780  IF  Ilav*1  THEN  Tsar- < T ida-Tv )» . 5 

2785  IF  I lav -2  THEN  Tsar-Tv 

2790  IF  Ihm-1  THEN 

2795  Tavg-7miet 

2300  Graa-37.9853+. 104988-Tavg 

2905  Taroo-HBS<Emt<7))*1 .2+6/ ( 1 0-Grad) 

2810  Tavgc*T in let-Tdrop- . 5 

2315  IF  ABS <  T avg- T avgc ) > . 0 1  THEN 

2820  Tavg- < Tavg+Tavgc > * . 5 

2325  GOTO  2800 

2830  END  IF 

'335 ' 

2340 1  COMPUTE  HATER  PROPERTIES 
2345  Kui-FNKw (  Tavg ) 

2350  Muua-- NMuiiH  Tavg) 

2855  Coui-^NCpijt  Tavg  > 

2360  Prij*- MP r'j(  Tavg ) 

2365  Rho(i/-i-NRhoij<  Tavg) 

2370  Uii-Tavg 
2375! 

2390!  Compute  MDOT 

2335  Mdo i-3 .9657E-3 +Fms— < 3.61 355E-3-Fms—<  3 . 820  OSE'G'Fms-* 1 . 23688E~ 7 -Fms-4 . 31397 
E- *  0  >  >  > 

2990!  **00 t  -Mdo  t  •<  1  .  0365-Tmiet-<1  .  966 44E-3-7 inlet *5.25 2E■6>>/,  .0037 
'’995  Mat*1) 

2900  Q-Mdo  t-Cpiii-Tgroo 

2905  Lntd* 'drop /LOG  ( <  T  mlet-Tsat)/<  T i  n i e  t- Tdrop-'sat) > 

29,0  Uo-i}/<PI-Oo*L-Lmca> 

29 1 5  Rw-Do-'.GG <  Do/  D  i )  /  <  2  .  cu) 

2920  Ty* T^at ++r-Ln td 

2925  Uio-^do r/ < Rhou-P I *0 1  2/4) 

2930!  IF  < c t * 0  THEN 
2935!  Tat- : 

2941)!  "droo- "drop*  .  004-Vij  2 
■xi/lC  i  ’!]+^  1 1  •  ■) 

2950!  2ND  "  IF 

2955  R»ii)-Ahoiu-')'j»0  ’.  /  Muu»a 

1260  4 1  ij/0  i  .  3 -A  T‘j  <'./?..>*<  Niidta/  F'lMu.u<  Fiji  )  )  .  '4 

2963  i c*  "avg— 5/  <  ?  I  -  0o-‘_ **i .  ’ 


-"6S;  j  i  -  ’  j  10  )  "•  .  J  ’  >.E:i 
2975  3w :  *  <  "y :  *  *y  i  .*  >  *  .  5 
_930  jQT.j  _9b0 
3935  END  IF 
2990  ~iii  i *  <  T'ji  *  ru  ic  >  •*  .  5 
2995  Ho* 1  /  ( 1  /Uo-Do/  < Q  i  -*H  i  >  -Ru  > 

2000  rhe  tac-OZ <  Ho-**5 1  •Do*'-  > 

2005  Ty * Tsat +*he tac 
20' 0  'he tao* "w- "sac 
3015  IF  Thetao<0  THEN 
3020  BEEP 

3025  INPUT  "TUALLCTSAT  (Q-COMTINUE.  !-£ND>'\Iev 

3030  IF  lev»0  THEN  GOTO  '325 

3035  IF  lev  I  THEM  3535 

3040  END  IF 

3045  END  IF 

3050' 

3055!  COMPUTE  VARIOUS  PROPERTIES 
3060  Tf  i  lm* <  Ty-Tsat  >  *  .  5 
3065  Rho*PNRho< r r i im ) 

3070  Mu*FNMu<  Tf  x In) 

3075  K-FNK<  Tf i im> 

3080  Co*PNCp<  Tf  i  lm ) 

3085  Beta-PNBeta( Tf i 1ml 
3090  Hf  g-PNHf  g<  Tgat ) 

3095  Ni*Hu-/Rho 
3100  Alpna-</<Rho*Cp> 

2105  Pr*Ni/Alpna 
3110  Psat-FNPsacI  Tsat ) 

3!  15! 

3120'  COMPUTE  NATURAL -CONVECTIVE  HEAT- TRANSFER  COEFFICIENT 
3125'  POR  UNENHANCED  END<S> 

3130  Hbar* ' 90 

3135  Fe«(Hbar-*,:,/<Kca*A)  >  ‘  ,5*Lu 

3140  Tann-PN Tanh ( Fe > 

3145  Theta* ‘hecao-'ann/Fe 

3150  Xx*<9 . 3 1 »9eta*ThetaD*0o  3*Tann/ < Fe*N i *Alpna> )  .166667 

3 ' 55  ry-< 1 *< ,559/Pr>  ( 9/ ’ 6 > >  <8/27 ) 

3'60  Hbarc-*'/Do'*<  .o*  .  387*<x/Ty  )  2 

3' S5  IF  ABS< < Hbar-Hbard /Hbarc ) > . 00 1  THEN 

3170  Hbar*<  Hbaf-Hbarc  >  •  .  5 

3175  GOTO  3135 

3190  END  IF 

3185 ! 

3190'  COMPUTE  HEAT  LOSS  RATE  THRCUGN  UNENHANCED  END<S> 

2'95  01  * ( Hbar*A*<cu»A >  . 5* The tac* Fann 

2200  0c*0-2*01 

3205  As-PI*D2*L 

?  ■?  <  Q  I 

32*5!  COMPUTE  ACTUAL  HEAT  FLUX  AND  BOILING  COEFFICIENT 

3220  Qdo-Gc/As 

2225  Htuxie-Gap/ 'hetaD 

1230  Cat  *<  Cp-*  The  tao/HTp  >/(  Qdp/ <  o  >*•<.  0 1  4/ <  9 . 3 1  *Hho  >  .5)  < 

''■'ft 

RECORD  '  I  HE  3F  DA'A  TAKING 
3245  Ic  Im-0  rHEN 
3250  OUTPUT  *09 : "TD” 


,3DD.2X.""3uik  Oil  X 


.  3D  .  3 . 5X . 


*  ih'r  *:  ^ in-.3 

-  "C  PQ  '*]  -T-  -  ' 

■;?o 

PRIM" 


2290 

CD 

►  R 

CO 

t— 

Cs 

n 

"OX.  "''Data  Set  Numoer  =  .  JDU  .  2X  .  "”’3uU  Ui- 

7,  -  1 

"".31 

3 . 2) .  5  X 

4A":.v 

.  3od  .  “o id's 

2295 

I  r  «.hw*0  HE) 

I 

2200  ^ 

PRINT  JS IMG  ' 

"0X.""’C  No:  1  2  2  -j  5 

2305 

3RINr  jSIMG  ' 

"0X.""'emp  3  (  tX.MDD.  3D  >";T<  J > . T<2> . T<3> . T 

(4)  . 

T<5) 

.  ’ «  S )  . 

23  i  0 

PRINT  'JSIMG 

"OX.’’”  'ja  Tlicd  'liad2  Tvapr  Psat  Tsujnp 

3315 

PRINT  JSING 

"1  OX  .2<MDD . DD . 1 X  > , ! X .HDD .  3D  .  1 X  ,2( 1 X  .HDD . 3D  > . 2X 

.  HDD 

.0"  : 

i  <ja «  >  i  ( 

r :  d2 . 

Tv  .Ps at .  Tslutip 

3320 

PRINT  USING 

'MUX.'"'  rhetaD  Htuce  Odp 

3325 

PRINT  USING 

'M  OX  .MOD.  3D  .  1  X  ,HZ.  3DE  .  !X  .MZ.3DE'' :  Thetao  .Htune  . 

Qop 

2330 

ELSE 

^  ^  cr 
JJJJ 

PRINT  USING 

'M0X.'"‘  Fms  Uij  Tsat  Tirti  Tdroo  ThetaD 

q 

Uo 

Ho 

3340* 

PRINT  USING 

'"OX ,4<2D.DD, !X> . Z. 3D . 'X .3D.DD. IX . 2<  HZ . 3DE . !X> 

’  ;  Fms  .  Vui 

.  Tsat . 

nlet . 

Taroo . Ihetac . 

Jo 

3345 

END  IF 

3350 

END  IF 

3355 

IF  lov*  t  THEM 

3360 

IF  J-’  THEN 

3365 

PRINT 

2370 

IF  Ihm-0  rh£ 

:  i 

3375 

PR INT  USING 

"  1  0 X  .  RUN  No  OiiOC  Tsat  Htooe  Odp 

The tao 

3360 

ELSE 

'HE 

2385 

PRINT  USING 

■••ox.""  F! is  oil;:  tsat  htube  odp 

~AB . 

3390 

END  IF 

3325 

END  IF 

3400 

IF  Ihm-0  'HEM 

3405 

PRINT  USING  " 

1 2X . 3D . 4X . DD . 2X . HDD . 3D . 3(  1 X . HZ . 3DE) " ; J . 3op . Tsat .H tuce . 3dP 

raa 

2J10 

ELSE 

24 ' 5 

=9 INT  JSIMG  " 

'  2X  .  3D  .  4X . 3D . 2X .HDD . DD . 3< : X .HZ .3DE>" :Fns . 3oo . 'sat .Htuoe . O' 

■w'J 

2420 

END  IF 

^42c 

END  IF 

24  30 

IF  In-0  'HEM 

3435 

3EE3 

2440 

INPUT  "OK  tq 

STORE  THIS  DATA  SET  < 1 -Y . 0-N) 7" . Ok 

3445 

END  IF 

3450 

IF  3k-:  3R  In 

*>  'HEN  J-J4-1 

3455 

IF  Ok  *  ’  PND  I 

n-0  'HEN 

3  46  0 

IF  Ihm-0  'HEN 

OUTPUT  »F • le i : Boo . To ioS . Ent  <  ♦ > . Vr . I r 

3465 

IF  Ihm*1  THEN 

OUTPUT  ?F  i  ie1  :6oo  .  ToidS.Emt  (  •  >  .Fns 

34  70 

ENO  IF 

34'S 

IF  la.:  -  '  *n£:l 

OUTPUT  iljf  .  ie  :  U.j  .  Jo 

3420 

IF  Ire-1  'HEN 

OUTPUT  aR  or  ^  1 0  ;  -  ns  .  RenJ 

3485 

IF  <  In  -  '  ]R  3 

k  »  ’  ;  JND  Ip  iot-'  'HEN  OUTPUT  i>P  1 0 1  :  OdD  .  The  tat) 

24'jO 

I  -  In-1  'HE'. 

3495 

BEE3 

25, ,0 

IMP'.!"  "NIL.  ■ 

HERE  BE  4N0THER  RUN  (  '  -  ' .  )-►!>  .  Go_nn 

l~  ■„ E 

Nr 'ir »  . 

Oo_ >n«0  'HEM  3525 

~  H  E !  I  ^  eo»  a  t 


f 


vjyl 

i 

I'A'I 

*&:{ 

?A 


m 

m 

I 

tyV 

i 

»v; 

Wi!| 

f 

lib 

;$S 

*5 

% 

v*» 


I 


•ji'i 

t 

i 


#!• 


f 


?5->;  SEE¬ 
DS  45  3RIMr  .'SING  "  '  MX  .  ""‘lO’E :  rata 

•  ,32_-‘  ;  i-5-’S 

3550  ASSIGN  »Fiiei  *0  * 

3555  GU'°Vr  *F  :  ie2  : Nrun- ' 

3560  ASSIGN  ?Filei  TQ  D)_*iie5 
3565  ENTER  “P  i  i  e  i  :  Date'5  .  ._d  tc  ’  .-jtc-.Itr 
35’0  OUTPUT  »F  :  ie2  :  Darei  .  Late  i  .  LatcC . 1 1 1 
3575  FOR  I*1  7Q  Nrun-' 

3530  IF  Ihtn-0  THEM 

3585  EMTER  *F :  le I : 3oo . ToioS . Eat ( * > . Vr . Ir 
3590  OUTPUT  $F  i  ie2  :  3op  .  To  ic)3  .  Emf  <  » )  .  Or  .  Ir 
3595  ELSE 

3600  ENTER  ?F  i  lei : 3oo . ToiaS .Erat <*>, Fms 
3605  OUTPUT  4F  i  ie2 : 3oo . To IdS . Emf < » > . Eras 
36  I  0  END  IF 
3615  NEXT  I 

3620  ASSIGN  »Fi  iel  TO  * 

3625  PURGE  "DUNMV“ 

3630  END  IF 
3635  BEE3 
3640  DR INT 

3645*  IF  Is  lot- 1  THEM  PRINT  USING  ”!  9X  .  ""NOTE : 

oiot  data  file  '”',1 
3650  ASSIGN  iFiie2  TQ  * 

3655  ASSIGN  iPLot  TO  - 

3660  IF  Iaf*i  rHEM  ASSIGN  ><J file  TO  * 

3665  IF  Ire*1  THEN  ASSIGN  iRe'iie  *0  * 


rata  runs  ye re  stored 


X-Y  pairs  yere  stored  in 


3670 

CALL  Stats 

3675 

3EE3 

2630 

INPUT  "LIKE 

rG  PLOT  DATA  ( 1 *Y , j *i 

3635 

IF  3k*1  THEN 

CAL.  3Iot 

3630 

SUBEMD 

3635’ 

3700! 

CURVE  ”ITS  OF  PROPERTY  -UNCTICMS 

3705 

DEr  -  NKcu<  '  > 

37'0! 

]F'-c  •:gp°cr 

250  '0  200  < 

3  7*5 

Tic  *  *  *273  .  *  5 

i  r  , 

2  720 

<*e?4-  .  '  T  2 - 

k 

37'": 

RE'URN  v 

3  730 

'NEMO 

-)  7  it; 

y  i  j 

DEr  PNt!u(T) 

3  74iV 

' 70  '3  360  < 

3'JRVE  'IT  OF  vrSCSU: 

3745 

■|c-r-Z"’3.  '5 

•  <:  *n  * 

3  750 

Mu-EXP < -4  .  j6 

26  *<  1  0  !  !  .  4  /  /  Tic )  >  *  ’  .01 

3755 

RE'URfl  Hu 

3760 

-MEMO 

3765 

DEr  rNCo  <  T > 

37 -T 

•30  '0  J00  < 

CURVE  rIT  OF  Ca 

"  7  7  ^ 

•»  r * _  •  c 

if  ~r  f 

Z?V) 

Ca-  •  4i) :  33*’  . 

6500  7E-  3*  '<  * '  .si4-4uE 

■  ■*OC 

Ia*0a*  '  )00 

1  'TQ 

PE 'URN  Ip 

i  7  ?  S 

-NEMO 

:3'>rj 

DEr  - MR6o 1 T  ) 

pone 

~v *'*;’ : .  -5 

if  -f  ,■ 

T  O  •  1 

>  '  .  •*"<.' 

753 .25  )  ,v  '0  R 

‘  •>  •  C 

3  o  *  ,5.  32'-)!  . 

’  4  5  4  ‘  4  *  ^  C  34  ’■1P>.*17 

J  '  I>  V_ 

~i  n  *  ►* » 3  /  . 

:2c 

-  E  JnN 

2-6 . 67353E-  1  0*T'<  3 


3341)  3  ICa<  7> 'EN«<  ~ ) 

3345  SE'iiRr!  =>- 
3350  "MEND 
3955  DEF  F  NK  <  T ) 

3360'  T <360  <  RITH  T  IN  C 
3365  K-  .  071 - . 00026 1  * 7 
3370  RETURN  < 

3875  EMEND 

3380  DEF  -N  T ann<  X ) 

3385  P-EXP(X) 

3890  O-t/P 

3395  Tann-<P-Q)/(P*Q) 

3900  RETURN  Tann 

3905  FNEND 

3910  DEF  FNTvsv(O) 

3915  COM  /Cc/  C(7) . Icai 
3320  T-C(O) 

3925  FOR  I-i  TO  7 
3930  T - E  *0 ( I ) *0 ' I 
3935  NEXT  I 
3940  IF  leal- 1  THEN 

3945  r-T-5.7422934E-2*r*<9.0277Q43E-3-T*<-9.32S9917E-5>  > 

3950  ELSE 

3955  T-T-8 .S26897E-2*T»<  3 . 76 ! 99E-3-T-5 . 0639259E-5) 

3960  END  IF 
3965  RErURN  T 
3970  FNEND 
3375  OEF  FNBecaO 
3900  RoD-^NRho ( r *  .  i  ) 

3985  Rom--NRho< T- . 1  ) 

3990  Beta--2/<Rop*Rom)*<Roo-Rom)/ .2 

3995  RETURN  3eta 

4000  FNEND 

4005  DEF  FNHtg<T> 

4010  Hfg-I , 3741 344E-5-E.< 3. 309436  IE -2*E-t .21651431 

40 ’5  RETURN  Hfg 

4020  FNEND 

4025  DEF  “MP-sat  (  Tg  ) 

40301  0  TO  30  3ea  -  CURVE  -IT  OF  3sac 
4035  Tf-i.3*Ec*32 

4040  Pa-5 . 345525t T t »<  .  : 5352082 <-Ef»<  ;  „4s3«0963E-3  +  rf*>}.6i5067!E-o> ) 
4045  Pg-3a- '4.7 

4050  IF  Pg>0  THEN  '  ♦•>38IG.-*;n  Hg 

4055  3aat-;3g 
4050  ELSE 

4065  3-3at-3g-29.32/'4.7 

4070  END  IF 
4075  PE'URM  3?ac 
4030  CNE'ID 

41185  DE-  -  NHsmoo  t  n  <  < .  3oo  .  I  sa;  > 

4090  DIM  a<5>  .3(5)  ..:<5,  .3(5  , 

4095  3RE4  . 20525 , . 25322 ..31 3048 . . 55322 . . 79909  .  '  .  00253 
4  1  00  DP  Er  .  745  •  5  .  .  72932  .  .  73' 39  .  .  ."225  .  .534  ’2  .  .  54  to- 
4105  DHEO  .4:092. 7726. .25142. .54H06. ,3'9'8. : .0345 
4  0  DPEA  .  7:  40  0 ,  .  / 231  j  .  .  725r>5  .  .  5 36o  ?1  ,  .  o65367  ,  .  6 1 339 

4 ' '5  READ  4<->.3<*>.5(»>.3<-> 

4' 20  .r  3oo'  -1  '4p‘l  .  -  o  o 

4 '25  IF  3oo *6  'HEN  I  ’4 


-i  ■  41.1  <P<  A<  ‘  )  «-0<  '  )  <  )  > 

-1 ■  :.Jc 

-i  ’  5  0  p  ( C  c  ’ )  *  Q  <  I  >  •*'_  j G  <  X )  > 

J'SE  END  :r 
Jib.)  RE'URN  -s 
J'bb  FNEND 


SUB  - 1 som  Of  .  C .. > 

COM  /W  i ;  /  D2 . D i . Do .  t.U.Hcu 
DIM  E.nf  <  : 2 ) 

ALISON  P[_GT  SUBROUTINE  DE 'ERMINE  OF  AND  Cl 
3EER 

INPUT  "ENTER  DATA  FILE  NAME" . Fi Le5 


INPUT  "ENTER  DATA 
BEEP 

PRINTER  IS  : 

PRINT  USING  "4X." 
PRINT  USING  "4X." 
PRINT  USING  "4X." 
PRINT  USING  "4X . " 
PRINT  USING  “4X . " 


Select  ootion: 

0  Vary  Cf  ana  Ci' 

1  z:x  Cf  Vary  Ci' 

2  Vary  Cf  Fix  C:' 

3  Fix  Cf  Fix  Ci1 


INPUT  "ENTER  0PTI0N".Icf ix 
PRINTER  IS  70! 

IF  Icfix-0  THEN  G535 

IF  Icf i x>0  THEN  BEEP 

IF  Icf  ix-1  THEM  INPUT  "ENTER  Cf.Csf 

IF  Icfix-2  THEM  INPUT  "ENTER  CI".Ci 

IF  Icfix-3  THEN  INPUT  "ENTER  Cf.  Ci" 


IF  Icfix-2  THEM  INPUT  "ENTER  CI".Ci 
IF  Icfix-3  THEN  INPUT  "ENTER  Cf.  Ci".Csf.Ci 
PRINTER  IS  ! 

INPUT  "Nant  To  Varv  Coef ??(*•*. 0«N)" . Iccoet 
IF  Iccoef-i  THEN  INPUT  "ENTER  COEr F".R 
PRINTER  IS  701 

IF  Icfix-0  OR  Icfix-2  THEN  Cfa-.004 
IF  Icf  ix*’  THEN  Cfa-Csf 
Ci a- Ci 
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APPENDIX  B 

FLOWMETER  CALIBRATION 

A  Fischer  Porter  Flowmeter  was  used  in  the  experimental 
apparatus  to  indicate  the  water  flowrate.  Prior  to 
conducting  the  boiling  tube  data  runs,  calibration  of  the 
flowmeter  was  performed  over  a  temperature  range  from  19  °C 
to  38  UC.  The  goals  of  the  calibration  were: 

1)  to  develop  a  mathematical  expression  relating  flowm¬ 
eter  percent  reading  to  mass  flowrate  and 

2)  to  determine  a  viscosity  correction  factor  due  to 
varying  temperature  range  and  flowrates. 

A  weigh  tank,  platform- type  scale,  and  stopwatch  were  used 
to  record  water  collection  data. 

Data  were  processed  using  the  program  "FMCAL",  listed 
later  in  this  appendix.  This  program  took  the  flowmeter 
percent  reading,  water  weight  collected,  and  elapsed  water 
collection  time  as  the  inputs.  The  outputs  were  an  experi¬ 
mental  mass  flowrate  and  a  mass  flowrate  difference.  This 
difference  was  given  by  a  mass  flowrate  computed  from  a 
correlation  minus  the  experimental  mass  flow  rate  for  the 
same  flowmeter  percent  reading.  The  correlation  was  based  on 
a  fourth-order  least  -  squares  fit  to  the  calibration  data. 

After  all  the  data  runs  were  completed,  the  data  points  were 
combined  into  one  file  and  a  single  fourth-order  polynomial 
was  generated  to  describe  the  relationship  between  the 
flowmeter  percent  reading  and  the  mass  flowrate.  Of  the 
original  data  points,  98.61  percent  were  within  a  ±5.5 
percent  range.  Reprocessing  the  data  points  led  to  78 
|  percent  within  a  ±  1.02  percent  range.  The  equation  (B.  1), 

resulting  from  the  calibration  data  analysis,  is  a  fourth- 
order  polynomial  from  a  curve  fit  of  the  reprocessed  data 
without  a  viscosity  correction  factor: 
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Experimental  discharge  rates  are  less  than  the  rates 
calculated  from  theoretical  flow  equations  primarily  because 
of  the  effect  of  internal  molecular  friction  of  the  fluid 
and  the  viscosity.  The  theoretical  flow  equation  is: 

r:SV-C3.-o  )  I  0  .  5 

’  -w-d  (B-2) 

where  Q=volumetric  flowrate,  A  = cross  -  sectional  area  of  the 
narrowest  part  of  annulus,  V.-svolume  of  float,  Pc=  density 
of  float,  Ac=cross-sectional  area  of  the  largest  part  of 
float,  and  c  =density  of  fluid.  The  Coefficient  of 
Discharge,  C^,  is  a  function  of  viscosity.  A  constant  value 
for  C.  in  equation  (3.2)  would  be  desirable  as  this  would 
show  a  negation  in  the  variations  due  to  the  viscosity.  The 
effect  of  the  float  shape  on  C .  with  varying  Reynolds  number 
is  shown  in  Figure  3.1  (a).  The  Reynolds  number  for  the 
flowmeter  uses  an  effective  diameter  of  the  tube  inner  diam¬ 
eter  minu3  tne  float  outer  diameter.  This  plumb  bob  style 
float  presents  a  large  surface  area  to  the  fluid  stream  and 
the  viscous  effects  increase  as  the  flow  rate  increases. 
[Ref.  31:  ?p.  9304-9803] 

The  data  collected  from  the  calibration  runs  showed 
little  or  no  change  in  the  Discharge  Coefficient  over  the  20 
to  90  percent  flowmeter  setting  range  and  the  15  K  change 
in  fluid  temperature.  The  square  edge  plumb  bob  float 
reaches  a  constant  Coefficient  of  Discharge  at  Reynolds 
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number  >  10000.  This  is  due  to  the  effective  shape  of  the 
float  on  the  streamline  pattern  as  shown  in  Figures  B. 2  (a) 
and  (b).  The  experimental  calibration  range  of  1300  <  Re  < 
27000  displayed  a  flattening  out  of  the  Coefficient  of 
Discharge.  A  different  effective  shape  is  seen  by  the 
streamline  pattern  because  of  the  minimal  clearances  between 
the  tube  bead  guides  and  the  float.  Also,  the  float  used  in 
this  particular  flowmeter  had  tapered  edges  versus  square 
edges.  This  different  effective  shape  is  shown  Figures  B. 2 
(c)  and  (d).  Figure  B. 1  (b)  shows  the  effect  of  the  float 
shape  on  the  Coefficient  of  Discharge  for  the  tapered  edge 
float. 

Additionally,  the  accuracy  of  the  flowmeter  contributes 
to  the  flattening  effect  on  the  Coefficient  of  Discharge 
shown  by  the  experimental  data.  The  accuracy  effect  is  felt 
in  two  major  components,  reproducibility  and  scale  factor. 
Reproducibility  accuracy  is  minimally  determined  by  the 
change  in  the  fluid  flowrate,  corresponding  float  movement, 
and  observation  of  float  position.  The  calibration  data  had 
good  reproducibility  over  the  various  flowrate  and  tempera¬ 
ture  ranges.  The  scale  factor  accuracy  is  inherent  to  the 
scale  markings  on  the  flow  tube.  Fischer  Porter  conducted 
numerous  tests  of  flowmeter  calibration  readings.  They  used 
a  maximum  scale  fraction  error  of  0. 9  mm  on  a  250  mm  tube. 
Based  on  their  results,  Fischer  Porter  reported  [Ref.  31: 
pp.  9814-9816]  an  accuracy  of  plus  or  minus  one  percent  at 
high  flow  rates  was  a  reasonable  expectation.  Also,  at  low 
flow  rates  ( <  20  percent)  a  resonable  expectation  would  be 
an  accuracy  of  eight  to  ten  perecent. 

The  experimental  data  followed  a  trend  of  negligible 
viscosity  effect.  The  flow  tube  design,  float  shape,  and 
scale  factor  accuracy  were  the  main  contributors  to  this 
negation  of  a  viscosity  correction  factor.  As  such,  equation 
(B.  1)  was  used  without  a  viscosity  correction  factor  to 
calculate  the  mass  flowrate  for  the  boiling  tube  runs. 
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APPENDIX  £ 
MODIFIED  WILSON  PLOT 


As  referenced  in  Chapters  III  and  IV,  this  appendix 
provides  the  data  reduction  procedures  used  to  compute 
inside  and  outside  heat-transfer  coefficients.  The  tubes 
used  for  data  collection  were  the  smooth,  modified  High 
Flux,  and  Modified  Turbo-3  tubes.  The  modified  Wilson  plot 
anaylsis  could  not  be  directly  run  on  the  enhanced  tubes 
since  reliable  boiling  correlations  for  these  surfaces  do 
not  exist.  Also,  only  two  unknowns  can  be  allowed  between 
correlations  for  the  inside  and  outside  heat- transfer  coef¬ 
ficients.  The  steps  used  in  this  procedeure  are  outlined 
below  and  are  contained  in  DRP4,  which  is  cited  in  Appendix 
A: 

Basic  Equations 


The 

overall 


overall 

thermal 


heat-transfer  coefficient 
resistance  is  given  by: 


in  terms 


of  the 


:j 

o 

where 


(C.  1) 


(C.  2) 


Modif led  Wilson  Plot 


1. 


Assume  C3 £  and  for  the  boiling  side  and  compute: 

rO.  S 

(C.  3) 


?rx  • 


,  (  o ,  -  o  )  i 

1  1  V 


U- 


_2 


i 

4 


h  = 


(csf) 


(C.4) 


Note^  with  r  equal^ to  0. 333 ,  equation  (C.3)  represents  the 
1  1  [Ref.  23:  p.  969]  . 


Rohsenow  correlation 
2.  Assume  C  and  compute: 
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5.  Compute  a  new  set  of  values  for  C.  and  C  -  as  follows: 

•  1  3: 


(C.  11) 


and 


C  :  C 
Sf 


(C.  12) 


6.  Repeat  steps  1  through  5  until  convergence  for  Cj,  and 
Cs/  between  two  successive  iterations  is  less  than  0.  1 
percent. 


APPENDIX  B 
DATA  RUNS 

The  following  table  outlines  the  data  runs  used  in  this 
investigation: 

Where: 

1.  WH  -  Smooth  tube 

2.  HF-  High  Flux  ( i. e.  ,  porous  -  coated  Korodense)  tube 

3.  HFM  -  Modified  High  Flux  ( i. e. ,  porous  coating  machined 
off)  tube 

4.  TB  -  Turbo-B  tube 

5.  TBM  -  Modified  Turbo-B  (i.e. ,  external  enhancement 
machined  off)  tube 

6.  T  -  inlet  water  temperature 

7.  V  -  water  velocity 

8.  const  -  constant 

9.  dec  -  decreasing 

10.  inc  -  increasing 


TABLE  4 

DATA  RUN  DESCRIPTION 
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U2IC 

ANALYSIS 

PZP.CENT3 

file 

WHO  5 

WHO  5 

HF13 

HF13 

T3C4 

T304 

,  q  . 

(wyn  ) 

( nx/  s) 

3 1140 

2.  24 

26690 

1.  39 

49540 

2.  24 

33730 

1.  39 

63  730 

1.  37 

53670 

1.  07 

5* 

* 

1.  56 

2.  51 

1.  5  6 

2.  51 

1.  43 

2.  51 

JRc 

Re 

1.  79 

2.  65 

1.  73 

2.  61 

1.  61 

2.  61 

5 a 

R 

1.  59  . 

2.  52 

1.  59 

2.  52 

1.  46 

2.  52 

Ji?m 

urrc 

4.  96 

3.  59 

3.  10 

2.  45 

2.  29 

1.  56 

JR 

4 

l 

2.  76 

2.  76 

11.  17 

11.  17 

2.  92 

2.  92 

3 

u 

5.  21 

4.39 

3.43 

3. 52 

2.  71 

3.  02 

3 

5h  . 

H. 

! 

2.  03 

2.  56 

1.  61 

2.  25 

1.  43 

2.  21 

5.  34 

5.  Q3 

4.  90 

4.  89 

3.  35 

3.  96 

Noce:  Ail  runs  were  perroraed  in  pure  rexrrgerant 
WH  ~  Smooch  rube 

Hr  =  Koroaense  cube  wich  High  Flux  Coating 
13  =  Turba-3  cube 

The  large  uncertainity  for  Che  wall  -resistance  of  file 
HF13  (?orou3-  coated  Korodense  Cube)  is  due  to  Che  10 
percent  uncertainly  of  Che  themal  conductivity  coeffi¬ 
cient.  For  this  investigation,  the  coefficient  used  was  an 
average  value  for  Che  values  given  in  Che  open  literature 
for  copper-nickel  places  or  Cubes. 
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APPENDIX  £ 

LIST  OF  NOMENCLATURE 

I.  NOMENCLATURE 

A  -  Surface  Area 

c_  -  Specific  heat  at  constant  pressure 

r . 

Cc  -  Oil  concentration 

-  Inside  Sieder-Tate- type  coefficient 

-  Rohsenow  coefficient 

d^  -  Inside  diameter  of  tube 

d0  -  Diameter  of  mouth  of  reentrant  cavity 

D  -  Diameter 

e  -  Depth  of  internal  ridge  or  fin 

Ev  -  Energy  transfer  per  unit  volume 

f  -  Friction  factor 

g  -  Acceleration  due  to  gravity 

g_  -  Gravitational  constant 

h  -  Co-nvective  heat  -  transfer  coefficient 

hc,-  Specific  enthalpy  of  vapor-liquid  mixture 

-  b 

k  -  Conductive  heat - trans f er  coefficient 

kT  -  Conductive  heat  -  transfer  coefficient  of  metal 

L  -  Tube  length 

L^  -  Entrance  effect  length 

.75  -  Mass  flow  rate 

p  -  Distance  between  ridge  or  fin  peaks 

Pr  -  Prandtl  number 

q  -  Heat  flux 

Q  -  Total  heat  transfer 

r  -  Tube  radius 

Re  -  Reynolds  number 

R..  -  Thermal  wall  resistance 

St  -  Stanton  number 
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U0  -  Overall  heat - t rans f er  coefficient 
Vj  -  Specific  volume  of  vapor 
V  -  Velocity 

A  -  Change  (  i.  e.,  Tin  -  Tout  ) 
o  -  Surface  tension 

-  Surface  tension  of  oil 

-  Surface  tension  of  refrigerant 
u  -  Viscosity 

uo  -  Viscosity  of  oil 

-  Viscosity  of  refrigerant 
p  -  Density 

oo  -  Density  of  oil 

Density  of  refrigerant 

2.  SUBSCRIPTS 
f  -  fluid 
i  -  inside 
in  -  inlet 
m  -  mixture 
o  -  outside 
ol  -  oil-liquid 
out-  outlet 

rl  -  refrigerant- liquid 

s  -  smooth  tube 

sat-  saturation 

w  -  water 

wo  -  outside  wall 
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ABSTRACT 


Heat- transfer  measurements  were  made  for  a  single, 
water-heated  tube  in  a  pool  of  R-114  to  simulate  operating 
conditions  of  a  water  chiller.  Data  were  obtained  for  a 
smooth  copper  tube,  and  for  two  commercially  available 
tubes:  a  spirally  roped  copper-nickel  tube  with  a  porous 
coating;  and  a  copper  tube  with  a  structured  outer  surface 
and  a  multiple-start  helical  ridged  inner  surface. 
Measurements  were  made  for  refrigerant-oil  mixtures  at  oil 
concentrations  from  0  to  6  mass  percent  with  a  boiling  pool 
temperature  of  13.  8  °C.  Results  for  the  two  enhanced  tubes 
with  and  without  oil  are  compared  to  the  smooth  tube  data. 
Enhancement  factors  for  the  overall  heat- transfer  coeffi¬ 
cient  were  4.0  and  3.6  for  the  structured  surface  and 
porous-coated  tubes,  respectively,  in  pure  refrigerant  and 
at  a  water  velocity  of  2  m/s.  For  these  same  conditions  the 
enhancement  factors  for  the  outside  heat- transfer  coeffi¬ 
cient  were  14.6  and  6.4  for  the  porous-coated  and  structured 
surface  tubes,  respectively. 


